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FOREWORD 

The meaSurement of atmospheric parameters to support activities meen 
on National Ranges involves the entire meteorological environment from Se 
the most Simple observation of data near the surface to the complex memes 
and difficult measurements of high altitude parameters by the use of sats 
rocket and satellite-borne instrumentation. The needs for accurate \ N 
assessment of the meteorological environment in which missile and 
Space system operations are conducted do not vary significantly from Ars 
one range to another, though the frequency and quantity of observa- ees 
tions may be considerably different. In recent years, many test Oyen 
programs involved the use of more than one range ard requived similar ste 


types of meteorological measurements from each range participating. 


.t 
*! 


4easurements of the atmospheric environment art, in most cases, 
ude by use of the same type of instruments on several ranges. To 
ensur2 the highest quality processed data, the Meteorology Group (MG) 
of the Range Commanders Council (RCC) is assigned the task of 
standardizing the methods used in the reduction of range meteorologi- 
tal data. A comnittee composed of the following members was formed to 
pursue 4a task titled "Standardization of Range Meteorological Data 


weety Tiare”; 
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1 H. c. Herring ESMC, Chairman 

Mr. F. J. Sehmidlin NASA/WFC 

tr. oa. Dunaway WSAR 

4r, £E. J. Keppel AD 

ar 3. a. BOlre WSMC 

ne J. A. Lea PMTC 

“Pie te E. Fasher AFSC 

fu, nn. 3. Noonan AFFTC 

" x. od, Evert NWC 

e x, Titus NTS/DOE 
“ 4. o£. Raffensburger UTTR 
Mr.oc, Glode@au KR 


sue Originas issee-ofPart I of this document was published as _ 
.t.a0 Gocumert 7 T08-67..in_ August 1967, Both Parts I and II of the 


} gen y were later published as IRIG Document 108-72. Tne current 
re,.s.on of the document is comprised of computer program documenta- 
Sacp for two different computer systems. The program in section A is 


uses by a Large central computer to process data from either 
nonrtranspoander or transponder rocketsondes. The program in section B 
is for the minicomputer (NOVA=-3/12) of the Meteorological Sounding 
osystea (MSS) und is used for processing data from the transponder 
rocketsondes. 


The U.S. standard computer program for the processing'of Robin 
Falling Sphere data to calculate winds, temperatures, pressures and 
densities is the 19/7 UDRI Robin Program. Other programs should be 
considered experimental and not used in lieu of the U.S. standard at 
the present time, In tne future, providing adequate comparison 
2xperiments are conducted on some of the newly developed falling 
sphere programs which may show improvements over the U.S. standard, 


the RCC MG committee on Standardization of Range Meteorological Data |. 


TN 


Reduction may recommend a new U.S. standard or approve other programs 
for use in addition to this standard for certain special applications. 


Range meteorological agencies are encouraged to conform to the 
standard data reduction techniques and to use the formulae contained 
in this document in an effort to achieve the highest possible degree 
of standardization of meteorological data at all ranges. In addition, 
the committee encourages other agencies involved in machine processing 
of meteorological data to conform to these standardized procedures to 
the maximum extent possible. The committee and the MG invite and 
solicit comments and recommendations for improvements in the data 
reduction techniques and associated information contained in this 
document. Pijeas2 forward such inputs to the Secretariat, Range 
Commanders Council, ATTN: STEWS-SA-R, White Sands Missiie Range, New 
Mexico 38002-5002. 
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: 
Be ABSTRACT 


- The Rocket Reduction Program is designed to process data acquired by standard Loki and Super 
Loki rocketsondes. This program accepts input data taken from the AN/TMQ-S5 record and 
the strip chart from the MSS Tracking System, a radar and/or AN/GMD-4 Rawin Set, plus co- 
rawinsonde baseline data. The program processes input data and outputs as its final product rocket- 
sonde data in readily usable tabular format; card decks, if desired; and rocketsonde message code 
(ROCOB). 
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10 MODULE DESCRIPTION 

1.1 Program ROCK 

Initializes data control arrays and controls flow of data through the subroutines. 
1.2 Subroutine SIMQ 


A least squares method of smoothing, provides a best fit of the components by generating a second~ 
degree polynomial equation. 


13 Subroutine OUTPUT 
‘Formats the data for output. 
1.4 Function WDIR 
Computes wind direction in polar meteorological coordinates and windspeed in meters/second. 
15 Function VALUE 
Sets missing data to nines or interpolates for value. 
1.6 Subroutine GRAVITY 


Computes the ratio of the station’s acceleration of gravity to gravity at 45 degrees latitude using 
the expression given at the top of Table 168, SMT, List 1968 (Bibliography Number 1). 


1.7 Subroutine CHKOUT 


Prints the tracking data in quality control format. Converts tracking data to Cartesian coordinates 
and writes data to a file. 


18 Subroutine RD46CH 
Reads the radar data and converts time to units and format consistent with TAER format. 
1.9 Subroutine TMPMERG 


Identifies type of temperature data, computes temperature from ordinate ratios, and combines 
temperature with tracking data versus time. 


1.10 Subroutine INTPTMP 


Identifies missing temperatures and interpolates to obtain temperatures through layers of less than 
3.0 km thickness. 


1.11 Subroutine MIDINT 


Computes the midpoint of altitude layers and calculates uncorrected component winds and fall rates 
for the layers between data points. Component winds and fall velocities are assigned to the midpoint 
of the layers. 
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1.12 Subroutine BANDTC 

Computes ballistic motion corrections for correcting winds according to Eddy (reference 3 in 
bibliography). Temperature corrections are computed for correcting temperatures as described in 
attachment 9 (reference 4 in bibliography). 

1.13 Subroutine COMPRSR 

Computes and tabulates the difference in temperatures from the co-rawinsonde observation and 
the rocketsonde observation in the overlap region. Chooses a baseline for barometric pressure 
computations. 

1.14 Subroutine TWOKM 

Smooths component wind layer data over 2-km thicknesses producing 2 km-mean layer winds. 


1.15 Subroutine PRESCOM 


Converts geometric altitude to geopotential units, computes mean virtual temperature, and computes 
the atmospheric pressure using a standard form of the hypsometric equation. 


1.16 Subroutine ONEKILO 

Rearranges data in a 1-km array and computes maximum altitude recorded by sounding. 
1.17 Subroutine PLOTCH 

Plots temperature versus altitude as an aid to quality control. 

1.18 Subroutine PUTOUT 

Outputs data to the printer file and the disk file. 

1.19 Subroutine ETFORM 


Formats data for entry into high altitude data formatting routine to develop data for archiving and 
publication. 


1.20 Subroutine CODIT 

Corrects units of temperature and converts other data. 

1.21. Subroutine RDRAWIN 

Reads co-rawinsonde data for input to the archiving and publication formatting program. 
1.22 Subroutine INT 


Interpolates for missing data in the input data. 
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1.23. Subroutine ROCOB 
— Encodes all rocketsonde parameters into the international ROCOB code and stores data for later 
| transmission. 


1.24 Subroutine MANDL 


Mandatory data level called by ROCOB to place data level into proper place in the code. 


i 20 


2.1 


where 


2.3.2 
where 
2.3.3 


where 


MATHEMATICAL DESCRIPTION 
Program ROCK 


Subroutine SIMQ 


Subroutine OUTPUT 


Computation of the Density of Each Level 


Q 
Q 
P 
T 


i 


348.38 - (P/T) 
density (grams/cubic meter) 
pressure of the level (millibars) 


temperature of the level (degrees K) 


Computation of the Velocity of Sound (Dry) 


331.45 + (1/273.15)'” 
velocity of sound (meters/second) 


temperature (degrees K) 


Computation of Vertical Wind Shear 


WS 


[((VX — VXP)? + (VY ~ VYP)’]'?/AALT 
vertical wind shear 

corrected east-west velocity component 

corrected north-south velocity component 
corrected previous east-west velocity component 
corrected previous north-south velocity component 


thickness of the shear layer 
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2.4 Function WDIR 
2.4.1 Computation of Wind Direction 
DC = tan” (VY/VX) 


If VX is positive, 


If VY is negative 


D= 
where D= 
DC = 
VX = 
VY = 


90 - DC 


270 -— DC 

wind direction, meteorological polar coordinates 
wind direction, Cartesian coordinates 

corrected east-west velocity component 


corrected north-south velocity component 


2.4.2 Computation of Windspeed 


C= 


where Cc 


VX 


VY 


(Vx? + vy)” 
windspeed (meters/second) 
corrected east-west velocity component 


corrected north-south velocity component 


25 Function VALUE 


None 


2.6 Subroutine GRAVITY 


2.6.1 Computation of the Ratio of Local Gravity to Gravity at 45 Degrees Latitude 


GRAT 


where GRAT 


$L 


1 — 0.0026373 - cos (2 - $1) + 0.0000059 - cos’ (2 - $1) 


ratio of local gravity to gravity at 45 degrees latitude 


station latitude (degrees) 


2.7 Subroutine CHKOUT 


2.7.1 rerhanling of Differences for Time, Azimuth, Elevation, Range, and Height Between 


Adjacent 
AG 


AZ 
Ag 


ve 
60/(GA — GB) 
(ZA — ZB) + AG 
(¢A — $B) - AG 
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a ete 
: ROCK 3 
a 
ms 
: A@ = (8A - 6B) - AG ch: 
Ry AR = (RA - RB) - AG 
i where AG = difference in time in fractional parts of minutes 
AZ = difference in geometric altitudes 
: 4¢ = difference in elevation angles 
oe A4@ = difference in azimuth angles bee 
! AR = difference in range ee 
GA = time of upper level (seconds) SS 
GB = time of lower level (seconds) mie 
ZA = altitude of upper level (meters) — 
ZB = altitude of lower level (meters) Re 
oA = elevation angle of upper level (degrees) am 
$B = elevation angle of lower level (degrees) Lae 
6A = azimuth angle of upper level (degrees) 3 oe 
6B = azimuth angle of lower level (degrees) Sy Ais 
RA = range of upper level (meters) ae 
RB = range of lower level (meters) ee 
2.7.2 Computation of Geometric Height << 
Z = ((RE + Hy)'+ R’ +2 - (RE + Hy) - sin $}'" — RE ee 
where Z = geometric height (meters) 
RE = radius of the Earth (meters) ee 
Hy = station height (meters) (see paragraph 7.4) en 
4 = elevation angle of the target (degrees) me 
R = slant range to target (meters) ma 
2.7.3. Correction of Elevation Angle for Curvature of the Earth re 
6. = 4 + {[Z + cos $/sin $]/ 2.2 - RE} - 
A-8 Be 


SPIEL SS 


where ¢, = corrected elevation angle (radians) 
¢ = elevation angle (radians) 
Z = geometric height (meters) 
RE = radius of the Earth (meters) 


2.7.4 Computation of the East-West and North-South Position Components for Each Data 


Point 
X= R -: cos $ - sin @ 
Y = R - cos $& + cos 6 
where X = .-west position component 
Y = north-south position component 
R = slant range to target (meters) 


$. = corrected elevation angle (degrees) 
@ = azimuth angle of the target (degrees) 
2.8 Subroutine RD46CH 
2.8.1 Computation of the Azimuth, Elevation, and Slant Range Factors 


9a = 0.125 - 6 + IRZ 
d¢ = 0.125 - $ + IRZ 
R’ = 8.0 - R + IRZ 
where 64 = azimuth factor 
¢¢ = elevation factor 
R’ = range factor 
9 = azimuth angle (degrees) 
¢ = elevation angle (degrees) 
R = range (meters) 
IRZ = timing factor 
2.8.2 Computation of Azimuth Angle, Elevation Angle, and Slant Range of the Data Point 
6 = 180 - Oy 


A-9 
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ROCK 3 
x 
: $= 180 - $ Ce 
R = 1.953125 - R” 
: where 6 = azimuth angle (degrees) 
¢ = elevation angle (degrees) 
: R = slant range (meters) 
a ¢¢ = elevation factor 
@, = azimuth factor 
R’ = range factor 
a 2.9 Subroutine TMPMERG 
R 2.9.1 Computation of Temperature (Steinhart Equation) 
a T = 1/[A + (B+ LR) +(C + LR) 
~ where §_T = temperature in degrees Kelvin 
LR = log of the thermister resistance value 
* A = ) calibration constants supplied by t 
g = } the instrument manufacturer 
: 2.9.2 Computation of Temperature (Layton-Clark Equation) 
- T = B/(LR - A)-C 
fa where T = temperature in degrees Kelvin 
LR = log of the thermister resistance value 
A= calibration constants supplied by 
: = } the instrument manufacturer 
2.10 Subroutine INTPTMP 
2.10.1 Computation of Interpolation Ratio 
K = (G,; — G)/(G, - G,) 
where K = interpolation ratio 
c G, = time of upper bounding level 
a 
mee 
A-10 ae 


oe. 
ra ways 
aida’. 


ee 


ane G, = time of lower bounding level 


‘ 
eV > bee T 
Pea ae 

Exe 
aX 


G = time of interpolated level 


2 2.10.2 Interpolation for Missing Temperatures 

= T=7T,-K-(-T) 

2 where T = interpolated temperature (degrees K) 

sf T, = temperature, upper bounding level (degrees K) 


T, = temperature, lower bounding level (degrees K) 


nw 
fi 


interpolation ratio 

2.11 Subroutine MIDINT 

2.11.1 Computation of the Midinterval VX and VY 
VX = (X - XA)/(GA — G) 
VY = (Y — YA)/(GA — G) 


where VX = midinterval east-west velocity component 
VY = midinterval north-south velocity component 
yy XA = upper bounding level east-west position component wine 
* X = lower bounding level east-west position component ce a 
. YA = upper bounding level north-south position component = 
. Y = lower bounding leve! north-south position component ee 
GA = time after launch to upper bounding level (seconds) we 
G = time after launch to lower bounding level (seconds) : ; 
2.11.2 Computation of Fall Rate as 
» VZ = (ZA - Z)/(GA - G) ite 
. where VZ = fall rate (km/seconds) eae 
st . aes 
o ZA = altitude, upper bounding level (km) Sonar 
a te ea 
: Z = altitude, lower bounding level (km) pele 
a GA = time after launch to upper bounding level (seconds) =e 
G = time after launch to lower bounding level (seconds) . Ee 
s ae 
< 2.11.3 Computation of Mean Temperature, Mean Height, and Mean Time Sey: 
a Ges TP = (T + TA)/2 Ye5 
= ee WANDS 
wie” 
eee 
B ray 
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ROCK 3 
H = (Z + ZA)/2 oa) 
GG = (G + GA)/2 
where TP = mean temperature 
TA = temperature, upper bounding level (degrees K) 
T = temperature, lower bounding level (degrees K) 
H = mean height 
ZA = height, upper bounding level (km) 
Z = height, lower bounding level (km) 
GG = mean time (seconds) 
GA = time, upper bounding level (seconds) 
G = time, lower bounding level (seconds) 
2.12 Subroutine BANDTC 
2.12.1 Computation of a Delta Time 
AG = G, - G, ee 


where AG = delta time 
G, = time of upper bounding level 
G, = time of lower bounding level 
2.12.2 Computation of Local Acceleration of Gravity 
GL = GRAT - 9.80616 
where GL 


local acceleration of gravity 


GRAT 


ratio of local gravity to gravity at 45 degrees latitude 
9.80616 = acceleration of gravity at 45 degrees latitude 

2.12.3 Computation of Ballistic Motion Corrections (reference 2, Bibliography) 

2.12.3.1 Correction of East-West and North-South Velocity Components 


VX, - [(VX3 — VX,)/AG] 
VX = VX, - ——— 
X= VX. ~ Tz, - VZ,)/4G)] + GL 


VY, = [(VY; = VY,)/AG] 


VY = VY: ~ (wz, - VZ,)/AG] + GL 
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where VX = corrected east-west velocity component aos 
VY = corrected north-south velocity component 2 
VX, = lower level east-west velocity component 
VX, = intermediate level east-west velocity component 
VX; = upper level east-west velocity component 
VY, = lower level north-south velocity component 
VY, = intermediate level north-south velocity component 
VY; = upper level north-south velocity component 
VZ, = fall rate, lower level 
VZ, = fall rate, upper level 
AG = delta time 
GL = local acceleration of gravity 5 ae 
2.12.3.2 Computation of Ventilation Velocity nae 
hte as 
V = [VZ,? + (VX - VX,)? + (VY - VY,)7]!” pens 
where V = ventilation velocity 
VZ, = intermediate level vertical velocity 
VX = corrected east-west velocity component 
VX, = intermediate level east-west velocity component 
VY = corrected north-south velocity component 
VY, = intermediate level north-south velocity component 


2.12.3.3 Computation of Aerodynamic Heating Correction 


where 


af Te las Soa) MSS 
eels 


AA 


Klz - V’ 
aerodynamic heating correction 


aerodynamic heating term vs altitude (see attachment 8, table 1) 


index into K1 table (20—70 km) 


ventilation velocity Te 
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2.12.3.4 Computation of Time Lag Correction ot 
AB = K2z - (T; — T,)/AG 

where AB = time lag correction 
K27 = time lag term vs altitude (see attachment 8, table 1) 


T,; = temperature of upper level (degrees K) 


ei i sa ee ee ee 


T, = temperature of lower level (degrees K) 
AG = delta time 
Z = index into the K2 table 
. 2.12.3.5 Computation of Radiation Heat Loss Correction 


AD = K3z - T* 
where AD = radiation heat loss correction 


K37 = radiation heat loss term vs altitude (see attachment 8, table 1) 


Z = index into the K3 table 
T = temperature of the thermistor (degrees K) . - 
2.12.3.6 Computation of Radiation and Electrical Correction 
AC = K4z 


AC = K57 if at night 


where AC 


radiation and electrical correction 


K47 = radiation and electrical loss term, day 


an 
a 

N 
{ 


= radiation and electrical loss term, night 


N 
Wt 


index into the K4 and KS tables 
2.12.3.7 Computation of Corrected Temperature . oe 
TC = T- AA + AB- AC + AD 


where TC 


corrected temperature 
T = temperature of the thermistor 
AA = aerodynamic heating correction 


AB = time lag correction 


sues 
| ong 
ae AC = radiation and electrical correction ae 
| ° AD = radiation heat loss correction bana 
2.13 Subroutine COMPRSR SS 
: 2.13.1 Computation of Interpolating Ratio Se 
K = (ZA ~ Zb)/(ZA — ZB) shan 
: where K = interpolating ratio a 
ZA = upper level height an : 
\ ZB = lower level height “ 
: Zb = baseline height oes 
. 2.13.2 Computation of Baseline Temperature meee 
: Tb = TA - [K « (TA — TB)] — 273.15 Bos 
: where Tb = baseline temperature (degrees C) oe 
: TA = upper level temperature (degrees K) 
; ee TB = lower level temperature (degrees K) of 
:: K = interpolation ratio : 
2.13.3 Computation of Temperature Difference 5 
; DIFF = Tb - TMP be 
: where DIFF = temperature difference (degrees C) . 
: Tb = baseline temperature (degrees C) : 
: TMP = co-rawinsonde temperature for the level (degrees C) Reve 
2.14 Subroutine TWOKM 
None sas 
2.15 Subroutine PRESCOM a 
: 2.15.1 Computation of Geopotential Height aa 
H = GRAT - [(RE - Z)/(RE + Z)] : aoe 
; where H = geopotential height bis 
GRAT = gravitational ratio, local gravity to gravity at 45 degrees latitude 
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RE = radius of the Earth Nee 
Z = geometric height 
2.15.2 Computation of Mean Virtual Temperature 


TV = (TV + TVP)/2 


l where TV = mean virtual temperature 


TV 


virtual temperature of the level 


TVP 


virtual temperature of the previous level 
2.15.3 Computation of Log of Pressure 
PL = PPL - (H - HP)/(PC - TV) 
where PL = log of pressure of the level 
: PPL = log of the baseline pressure 
PC = 67.442 = metric constant 
TV = mean virtual temperature of the layer 
i H = geopotential height of the level ‘. 
HP = geopotential height of the previous level 
2.16 Subroutine ONEKILO 
; 2.16.1 Computation of Next Output Level 


SY RAPE aoe 
ye 


, ee, # 


B ‘ef 


. 
«4, 


Hs = [(IH/1000) - 1000] + HINT 


where Hs = desired output level 


eens 
° 


IH = height converted to integer e 


HINT 


output interval requested 
. 2.16.2 Computation of Interpolation Ratio 
K 


(Hs - HB)/(HA — HB) 


where K interpolation ratio 
Hs = output level 


2 HA = height of upper bounding level 


HB = height of lower bounding level soe 
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2.16.3 Interpolation of Position Components, Temperature, Uncorrected Temperature, Fail Rate, 
and Log of Pressure 


Xs = K - (XA - XB) + XB rece 
where Xs = interpolated value wae 
K = interpolation ratio ie a 

XA = upper level value 

XB = lower level value 


2.17 Subroutine PLOTCH 

None 

2.18 Subroutine PUTOUT 

None 

2.19 Subroutine ETFORM 

2.19.1 Computation of Interpolation Ratio for Significant Data 


K = (Hs — HB)/(HA — HB) 


e where K = interpolation ratio a 
Hs = output level 

HA = height of upper bounding level marae 

eee 


HB = height of lower bounding level 


2.19.2 Computation of Interpolation Ratio for Mandatory Data ee 
K = (PL — PB)/(PA — PB) “ 
where K = interpolation ratio ~~ a 
PL = log of pressure at output level ee 

re 

PA = log of pressure, upper bounding level oe 

PB = log of pressure, lower bounding level hd 
ae 

2.19.3 Interpolation for Output Data o a 
5 Ae 

Xs = K - (XA - XB) + XB pee 

Aaa 

pe where Xs = interpolated value ~ oe, 


. . . » . . . > . - . . . te . ce - - SS te s = cee CA. ae Nec. ae : z Nee a - . . - # 
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K = 


interpolation ratio at, 
XA = upper level value 
XB 
2.20 Subroutine CODIT 
2.20.1 Computation of Velocity of Sound 


lower level value 


Vs = 331.45 + (1/273.15)” 
where Vs = velocity of sound (meters/second) 
T = temperature (degrees K) 


2.20.2 Computation of Pressure 

P = 10PL 
where P = barometric pressure (millibars) 

PL = log of pressure 

2.20.3 Computation of Density 

@ = (P/T) - 348.38 % : 

where @ = density (grams/cubic meter) 
P = pressure (millibars) 


T 


W 


temperature (degrees K) 
2.21. Subroutine RDRAWIN 
None 
2.22 Subroutine INT 
2.22.1 Interpolation for Output Data 
K = (Hs — HB)/(HA — HB) 
where K = interpolation ratio 
Hs = output level 
HA = height of upper bounding level 


HB = height of lower bounding Jevel 
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Xs = K - (XA — XB) + XB 


l | where Xs = interpolated value 
: K = interpolation ratio 
XA = upper level value 
i XB = lower level value 


2.23 Subroutine ROCOB 


None 

i 2.24 Subroutine MANDL 
None 
3.0 INPUT 


) Input to the rocket program is divided into four groups: temperature data, tracking data, co- 
rawinsonde data, and control data from the terminal. 


3.1 Tracking Data 


The program reads two basic formats of input data: the 46-character radar data and the TAER 
format. The other input formats are converted by editor programs to the TAER format. 


Rei 
ee 


3.1.1 Radar Tracking Data, 46-character 


The 46~character data is from the tracking radar and is put into a file in the following format 
(see attachment 6): 


7911214124802014250255460000 1 3634434622462 


See attachment 8 for a breakdown of this data. 


3.1.2. TAER Format and MATTY File are 


The TAER format and MATTY file are identical except for the first line of the MATTY file (see 
attachment §). 


GGGGG.G AAA.AAA EE.EEE RRRRR.R 

" GGGGG.G = time of data (seconds after midnight) . ® 
: AAA.AAA = azimuth angle (degrees from true north) 
EF.EEE = elevation angle (degrees from horizontal) 


RRRRR.R slant range (yards) 


3.2 Temperature Data 


Temperature data is in three types of record groups collected by one of two types of tracking systems. 
The tracking systems are the AN/GMD-4 and the MSS system. The three types of records follow. 


3.2.1 Ordinate Ratios (Attachment 2) 
Record Contents 


1 Test number (5 digits) 
Station number (2 digits) 
Month (2 digits) 
Day (2 digits) 
Year (2 digits) 
Time of launch, GMT (4 digits) 
Type of instrument (3 digits) 
Plain language name of motor 


Type 
RAEMG group (see attachment 1, tables | through 5) (6 digits) 
3 Temperature equation constants 
Three 8-character numbers 
4—24 Temperature sensor calibration data 
Two 8-character numbers each 
25 Separator flag 
26—end Ordinate ratio (3 digits) referenced to 95 followed by time in seconds after 


launch (4 digits) 
3.2.2 Coded Temperatures (Attachment 3) 


Record Contents 
1 Test number (5 digits) 
2 Station number (2 digits) 


Month (2 digits) 
Day (2 digits) 
Year (2 digits) 
Time of launch, GMT (4 digits) 
Type of instrument (3 digits) 
Plain language name of motor 
Type 
RAEMG group (see attachment 1, tables 1 through 5) (6 digits) 
3—end Coded temperatures (4 digits) followed by time in seconds after launch (4 
digits). Temperature is coded by adding 50 to all negative temperatures and 
omitting the minus sign. 


3.2.3 Missing Temperatures 


Records | and 2 are the same as in the preceding paragraphs; in this case, they constitute the entire file. 
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ee as 3.3 Co-Rawinsonde Data 


The co-rawinsonde input data is the 250-meter output data from the rawinsonde program. Attachment 
15 illustrates the format of this file. Its contents are as follows: 


ae Altitude Meters 
~ Wind direction Degrees 
y:, Wind speed Meters / second 
Temperature Degrees Celsius 
Dewpoint Degrees Celsius 
Pressure Millibars 
Relative humidity Percent 
Absolute humidity Grams/cubic meter 
- Density Grams/cubic meter 
: Refractive index Microwave in N units, optical fs 
ei Velocity of sound Knots ba 
rs Shear Seconds” ws 
2. Vapor pressure Millibars es 
: Precipitable water Millimeters o 
ne 3.4 Control Data (Operator Dialogue) ; 
= s A sample of the operator dialogue follows: 
ye IS THIS A NIGHT RUN ? YES/NO 
a ? NO 
~ IS THIS A MATTY RUN ? YES/NO 
ce ? YES 
4.0 OUTPUT 
Ree The primary output of the program is constant altitude data in kilometer intervals, 250-meter intervals, 
te and 1,000-ft intervals. This information is output to the printer, disc file (or magnetic tape), and 
oe, punch card file. The format for this data is the same to all three files as follows: 
Height Geometric feet or meters at standard intervals of 1 kilometer, 250 
meters, and 1,000 ft 
Wind direction Degrees from true north. Missing data is shown as 999 
Windspeed Knots or meters/second. Calm is 000 and missing is 999 
Temperature Degrees Kelvin 
z Temperature correction Degrees Kelvin 
7 Atmospheric pressure Millibars 
ie oe Atmospheric density Grams/cubic meter 
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Fall Rate Feet /second or meters/second (determined by height units) 
Speed of Sound Knots or meters/second (determined by height units) 
Wind Shear Seconds” 


Additionally, in the printer output, there are three quality control sets of data and a rocketsonde 
message code. The first quality control data set contains tracking information as follows: 


OBB ed 4 ABB et 


Data point time Seconds after midnight 

i Azimuth angle Degrees from true north 

ES Elevation angle Degrees from horizontal 

Range Meters 
Height Meters 
Data point time Time of day — hours, minutes, and seconds 
Time difference Time difference between adjacent data points, seconds 
Azimuth difference Change in azimuth angle/minute, degrees 
Elevation difference Change in elevation angle/minute, degrees 
Range difference Change in range/minute, meters 
Height difference Change in height/minute, meters 


The second quality control data set contains temperature information as follows: 


Height Meters 

Temperature Temperature of the rocketsonde instrument, degrees Celsius 
Temperature Temperature of the rawinsonde instrument, degrees Celsius 
Temperature difference Difference between above two temperatures at the height indicated 


The third quality control data set is a graph of temperature in degrees Kelvin versus height in meters. 
The graph covers a range of 20,000 to 70,000 meters in height. 


The teletype output is the rocketsonde coded message. Attachment 10 shows a sample copy of 
this output and attachment !1 gives a breakdown of the code. 


The High Altitude Meteorological Data (HAMDATA) output is saved on disk for use as input to 


the HAMD program. A sample copy of this data is shown in attachment 14. The first section 
of this file consists of: ; 


vy 

ete e, 
tf 
whee 


vr 
e 
eae 
4 afr, i 


ie 
Tt 
Pd 


DRT 


nf 
: Column 
\ Sar 
a 1—5 
if 6—8 
x 9-11 
ie 
. 12—15 
16—19 
Z 20—23 
24—27 
28—30 
c 31-34 
. 35—37 
r e 38—44 
45—S] 
52—54 


Altitude in geometric decameters 
Wind direction with respect to true north, degrees 


Windspeed in whole meters/second 


Uncorrected y component, with a minus sign indicates northerly component, 
blank for southerly component 


Uncorrected x component, with a minus sign indicates easterly component, 
blank for westerly component 


Corrected y component, with a minus sign indicates a northerly component, 
blank for southerly component 


Corrected x component, with a minus sign indicates an easterly component, 
blank for westerly component 


Fall rate in whole meters/second 
Temperature in whole degrees Celsius 


Temperature correction, with a minus sign indicates negative temperature 
correction, blank indicates positive correction, whole degrees 


Pressure in millibars to four significant digits with one digit to the left of 
the decimal point, three digits to the right of the decimal point, and a signed 
“‘exponent’’ to indicate the power of 10 by which the number is to be 
multiplied. 


Col. 38: first significant digit (never zero) 

Col. 39: decimal point 

Col. 40—42: three additional significant digits 

Col. 43: sign of exponent (never blank) 

Col. 44: power of 10 by which number is to be multiplied 


Density in grams/cubic meter, same format as pressure 


Speed of sound in whole meters/second 


The second section of this file contains data as follows: 


re Height Meters 
Wind direction Degrees from true north 
Wind speed Meters / second 
Temperature Degrees Celsius 
Dew point Degrees Celsius 
Barometric pressure Millibars 

eS hae Refractive index N units 
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5.0 ALGORITHMS Ce 
5.1 Subroutine MANDL 
Convert mandatory level pressure to whole number for output in rocketsonde code using 
ws = (10PLO + Fry) 10)S/10) 
where IPS = pressure at mandatory level encoded as a whole number 


PL(J) = log of the pressure at the mandatory level being converted to 
a whole number 


Fl = round off factor used to move decimal point, this factor is 
based on PL(J) 


If PL) < 0, Fl = 0.05 
< -1, Fl = 0.005 
< -2, Fl = 0.0005 
< —-3, Fl = 0.00005 
< —4, Fl = 0.000005 


JS = intermediate factor based on the pressure level indicator J 
which is set by the number of levels processed 


if J > 31, JS 
< 31, JS 


—NWAAD 


5.2 Subroutine TWOKM 


Set up values that are to be used when calling SIMQ by multiplying elements of the matrix and 
assigning to proper place in symmetric coefficient matrix (A). 


AX, J) = XAq, yy + XK, y+ SSK, Fy 
AY¢, 5) = YAa, x + YC, > YOu, yy 


Set up values that are to be used when calling SIMQ by multiplying elements of the matrix and 
assigning to proper place in constant vector matrix (B) using 


Using the coefficients returned from SIMQ, calculate a smoothed x and y from 
SX = a'XB, + aXB, + XB, 
SY = a'XB, + aYB, + YB; 
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where a=5 
XB = coefficients returned from SIMQ 
SX = smoothed east-west component 
SY = smoothed north-south component 
5.4 Subroutine SIMQ 


The least squares method provides a best fit of the components by generating a second degree 
polynomial equation, of general form: 


y = kx’ + kx + ko 


and computing an offset from the center of the 9-point range. Two matrices are created designated 
{A] and [B]. Symmetric coefficient matrix (SCM) [A] is in the form: 


N =X, =Xj 

[A] = =X; =X} EX} 

EX} EX} EX} 

and constant vector matrix (CVM) [(B] is in the form of: 
i 


[B} ag EX; Yj 


ut 


2 

=X; Y; 
Note: < indicates summation from 1 to N. N=9 
Basic matric equation [A] [K] = [B] is solved by finding the inverse of [A], [A]"' and performing 
multiplications: 

[AJ [A] [K] = [AJ [B} 

{1} (K} = [A] [B} 

[A] | is calculated by SIMQ using the GAUSS-JORDAN inplace matrix inversion method. SIMQ 
returns the appropriate coefficients k,, k,, and ko, which the calling routine incorporates in calculating 
the appropriate values. 


6.0 FLOWCHARTS 


See the following pages. 
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6.2 Subroutine SIMQ 


COMPUTE 
DATA AND 
SAVE 


HOUSEKEEPING 


MAKE 
DATA 
COMPARISON 


RESET 
VARIABLES 


ABS (TEMP) 
<0 
? 
NO 


RESET 
VARIABLES 


COMPUTE 
NEW 
DATA 


COMPUTE 
TEMP | 


COMPUTE 
NEW 
VARIABLES 


COMPUTE 
INDEX 
Ic 


RETURN 
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> 3.1 Subroutine OUTPUT 


7 : 
% 

Ls 

= fs 
) Neat 
5 

. 


TiNar 
ware 

‘ by 
wah 
Pa Are 


PRESSURE 
MISSING 


SET 
PRESSURE AND 


: COMPUTE DENSITY MISSING count s oe 
~ DIR AND PRESSURE eee 
: SPEED ee 
, eee 
i SET : 

i DENSITY AND 

. SPEED OF 


SOUND MISSING 


CONVERT 
SPEED TO 
PROPER UNITS 


PRESSURE 
MISSING 


COMPUTE 
DENSITY AND 
SPEED OF 
SOUND 


SET TEMP 
AND 
VIRTUAL TEMP 
MISSING 
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6.3.2 Subroutine OUTPUT aE Cs 


WRITE 
HEIGHT 
OFRECTION 
SPEED 


CHECK 
BASELINE 
DATA 


RESET 
VARIABLES 


RETURN 
"ae Nee 
| on Are 
a VARIABLES : 
HEIGHT 
< 
20000 ? 
HEIGHT 
< 
65616? 
; a 
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6&4“ Function WDIR 


ENTRY WSPD 


SEY WIND 
SPEED 
MISSING 


COMPUTE 
WIND DIRECTION 


oes SUBTRACT ares 
le DIRECTION COMPUTE sae 


SUBTRACT FROM 270 
DIRECTION 
FROM 90 


WIND SPEED 


Se ee 
DEC ee ee ae 


65 Function VALUE 


+s 


ales 


COMPUTE 
VALUE 
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6.6 


Subroutine GRAVITY 


INITIALIZE 
CONSTANTS 


CHECK STATION 
NUMBER 
1F CORRECT 


CORRECT NUMBER 


READ IN 
LATITUDE AND 
LONGITUDE 


COMPUTE 
GRAT 


RETURN 
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6.7 Subroutine CHKOUT Ae 
COMPUTE 
CORRECTED 
Sraee POSITION 
COORDINATES 
SEF COMPUTE 
CONSTANTS one 
[cas oare \ 7 
bi Q.C. DATA 
ks \ gate time _// 70 PRINTER 
WRITE 
WAITE TRACK DATA 
TITLE DATA sor aenated 
TO PRINTER aches 
READ 
FIRST TRACK 
DATA CARD 
RADAR ves SET 
DATA ITS38 = 1 
? 
YES 
e 
: ENDFILE REWIND 
SCRATCH SCRATCH 
. FILE FILES 
NO 
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6.9.1 Subroutine TMPMERG 
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6.9.2 Subroutine TMPMERG 
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6.14 Subroutine TWOKM 
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6.18 Subroutine PUTOUT 
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6.20 Subroutine CODIT 
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Subroutine RDRAWIN 
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6.22 Subroutine INT 
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6.24.1 Subroutine MANDL 
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6.24.2 Subroutine MANDL 
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MATHEMATICAL AND PHYSICAL CONSTANTS 


Program ROCK 


6371229.315 
67.442 


(Mean Average) radius of the Earth, meters 
Constant in the hypsometric equation, SMT, List 1968 


See attachment 9, table 1 for constants ATC, BTC, CTC, CTD, and DTC. 


See attachment 12 for logs of mandatory pressures. 


7.2 
None 
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7.4 


ie 


None 


7.6 


1d 


Subroutine SIMQ 


Subroutine OUTPUT 
348.38 


331.45 
273.15 


Function WDIR 
57.29578 
Function VALUE 


Subroutine GRAVITY 


5 

4 

79 

28.29 

17.7 

7.58 
57.29578 
0.0026373 
0.0000059 


Subroutine CHKOUT 


0.9144 
2:2 


Gas constant for dry air with pressure in millibars and density 
in grams/cubic meter 

Speed of sound at zero degree Celsius, miles/second 
Converts degrees Celsius to degrees Kelvin 


Degrees / radian 


Station height of Cape Canaveral AFS in meters 
Station height of Antigua AAFB in meters 
Station height of Ascension AAFB in meters 
Latitude of Cape Canaveral AFS, degrees 
Latitude of Antigua AAFB, degrees 

Latitude of Ascension AAFB, degrees (South) 
Degrees / radian 

Gravity ratio constants from SMT, equation 

for Table 168 (Bibliography No. 1) 


Yards to meters conversion factor 
Constant used in Earth’s curvature correction to elevation angle 
(Bibliography No. 2) 
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7.10 
None 
TAL 
None 


7.12 


7.13 
None 
7.14 
None 
715 
None 


7.16 


TAT 
None 
- 7.18 


None 


Subroutine RD46CH 
86400 

3600 

60 

0.125 

180 

8 

1.953125 

90 

Subroutine TMPMERG 


3600 
60 


Subroutine INTPTMP 


Subroutine MIDINT 


Subroutine BANDTC 


9.80616 


Subroutine COMPRSR 


Subroutine TWOKM 


Subroutine PRESCOM 


Subroutine ONEKILO 


1.94254 
3.280833 


Subroutine PLOTCH 


Subroutine PUTOUT 


Number of seconds in a day 
Number of seconds in an hour 
Number of seconds in a minute 


Decodes 46-character 
Radar encoded data 


Number of seconds in an hour sates 
Number of seconds in a minute a 


Acceleration of gravity at 45 degrees latitude 


Meters/second to knots conversion factor 
Number of feet in a meter 
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7.19 Subroutine ETFORM 


None 

7.20 Subroutine CODIT 
None 

7.21 Subroutine RDRAWIN 
None 

7.22 Subroutine INT 
None 

7.23 Subroutine ROCOB 
None 

7.24 Subroutine MANDL 
None 
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~ 
a RAEMG CODE Tables: 
Table 1 
. Symbol R = reason for no report and ground equipment employed 
re, Code 
. Figures 
0 Launch not scheduled 
1 Rocket motor failure 
2 Instrument (or telemetry) signal not received 
3 Ground tracking equipment failure 
4 Range restriction prohibited launch 
5 Weather prohibited launch 
6 Lack of expendables prohibited launch 
7 Radar only employed 
8 Radar and telemetry equipment employed 
9 Telemetry equipment employed 
e Tables 2 and 5 
Symbol A = thermodynamic correction technique 
or 
G = wind correction technique 
e Table 3 
Symbol E = method of reducing data 
Code 
Figures 
0 Manually — nomogram 
Ae I Electronic computer 
os 2—8 Unassigned 
é 9 Other method 
Table 4 
@ 
ae Symbol M = type of wind sensing equipment 
Code 
Figures 
e 00 No wind sensor 
01 Chaff, experimental 
02 Chaff, metalized 
- 03—09 Unassigned 
: 10 Parachute, experimental 
11 Parachute, 0.5 to 3.5 meters diameter 
e 12 Parachute, 3.6 to 5.5 meters diameter 
“s : 13 Parachute, greater than 5.5 meters diameter 
a 14 Mesh decelerator, experimental 
30 Starute, experimental 
31 Starute, 0.5 to 3.5 meters diameter 
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Table 4 (continued) 


Code 
Figures 


32 
33 
50 
51 
55 
60 
65 
70 


ROCK 3 


Starute, 3.6 to 5.5 meters diameter 
Starute, greater than 5.5 meters diameter 
Sphere, experimental 

Sphere, inflatable 

Grenade, experimental 

Chemical trail, experimental 

Meteor trail, experimental 

Remote sensing, experimental 
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ATTACHMENT 3 
SAMPLE COPY TEMPERATURE DATA WITH CODED TEMPERATURE VALUES 
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BREAKDOWN OF 46-CHARACTER DATA 
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DATA FORMAT Stes 

COMPUTER OCTAL SEMI-RAW TTY OUTPUT ae 

9.” 4, 

oats 

(38/46 CHARACTER) Peay 

ey 

nan 

SOM Line: JU} <4 or DD! <=} Description 

Char. Char. or : : 
No. Range Content Se 
Vehicle ID { 
1 Line Feed Beacon Skin ' 
2 ft Figure Shift : ar 
3 uv Figure Sh:‘t 0 & i 

+ 9) Vehicle 1D R$»; 4 eat 

s (0-9) Station ID ! : wes 

6 (0-9) Station 1D | Pre a 

7 @-9) Radar ID ras 

g (OQ or 2) Validity ID ran 

9 (0-2) Hrs (tens) Ni , 

10 (0-9) Hrs (units) eae 


tH (0-5) Min (tens) 

12 (0-9) Min (units) 
13 (0-5) Sec (tens) 

14 (0-9) Sec (units) 

18 (0-7) Az 21, 20, 19 
16 (0-7) Az 18, 17, 16 
\7 (0-7) Az 15, 14, 13 


18 (0-7) Az 12, II, 10 5 

19 (7) Az 9, 8,7 Station ID 

20 (0-7) Az 6, 5, 4 

ny AST oe 
Apu (0-7) El 21, 20, 19 CRO ee 3 
23 (0-7) EI 18, 17, 16 HAW -- 5 
24 (0-7) El 1S, 14, 13 PAT 0.18 3 
Pan! (0-7) Ei t2, 11, 10 GBI 3.13 7 
2% (7) E1987 GTK 7.18 4 
2? (0-7) EI 6, 5, 4 REDSTONE 5 
28 (0-7) El 3, 2, 1 MLA 19.18 4 
30 (0-7) Rng 24, 23, 22 ANT 91.18 3 
31 (0-7) Rng 21, 20, 19 CNV 1 16 , 
2 (0-7) Rng 18, 17, 16 . 

33 (0-7) Rng 15, 14, 13 

34 (0-7) Rng 12. 11, 10 

38 (0-7) Rng 9, 8, 7 

36 (0-7) Rng 6, §, 4 

x“ (0-7) Rng 3. 2, 1 

38 < Carriage return for end of 38 character data 


38 (0-3) 46 char data day of year (100) 

x9 (0-3) 46 char data day o. year (10) 

40 (0-9) 46 char data day of vear (1) 

41 (0-9) 46 char data catalog number (1000) 
42 (0-9) 46 char data catalog number (100) 


43 (0-9) 46 char data catalog number (i0) 

44 (0-9) 46 char data catalog number (1) 

4s (0-9) SPACE 

46 < Carriage return for end ot 46 character data 


EOM Line: 3 *%s ¢ 


The information contained in this publication ts proprietary to Pan Am World Services Inc. It mav not be copied or 
reproduced except for reproduction by the U.S. Government for L.S Government purposes without the express permis- 
sion of Pan Am. 
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ATTACHMENT 9 
DESCRIPTION OF TEMPERATURE CORRECTION TECHNIQUE 
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d ROCK 3 
cS 
a 
a a (This attachment was extracted from Federal Meteorological Handbook 10) 
Description of corrections applied to observed Rocketsonde Temperature 
Temperature Corrections for Datasonde Instrument 
on 
t The Datasonde instrument employes a 10-mil bead thermistor, mounted on a loop of silverized 
» mylar, commonly referred to as the loop mount. The general procedures for applying temperature 
corrections to film mounted thermistors, as described by Henry (Bibliography Ref. 3) have been 
expanded and adapted to the specific configuration of the loop mounted thermistor as used in the 
Datasonde instrument. This work has been done by Krumins and Lyons (Bibliography Ref. 4), 
and others, anc has been adopted as the U.S. standard correction technique. The corrections 
developed for application to the loop mounted thermistor include acceptable values for incoming 
radiation heating corrections. 
ae 
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The true atmospheric temperature, Tatm> iS equal to the observed thermistor temperature, T,, plus 


certain correction values described by the following terms: (some rearrangement of the elements 
of each term has been made from the similar terms used for the Arcasonde instrument and the 
elements for the incoming radiation corrections have been added to the fourth term.): 


Aerodynamic heating term, 


« 
A;hy + weAghe ac. 


Time lag term, 


WeCeMe | aT, 


wages tp Pete 
air a 2 


[ em, ‘i 
K = + | 
j Argh, + weAghe dy 
L 
Radiation heat loss, term, 
Ele + WeAce 
kK, = + | cet SrA] oe 
Ayhy + weAghe 
Incoming radiation and electric heating term, 
a Ky = - ___ 2st Gee J 
‘ Ah, + weAghe 
wraeeG a4G WraceG 
= f"sfVsf J_ stVst ig £"sfVsf Alb. J 
Ah, + weAghe Ayhy + weAghe Ayh, + weArhe 
a1Gyy + Weay¢Gi¢ 4 Wr 
_————— oTerp - eS 
Ah, + weAghe Ayht + wpAphe 


The four correction terms applied to the observed thermistor temperature can be simply identified 
in a correction expression as follows: 


dT 
Tatm = Ty ~ (Ky(2V’) + kno] + (K3(Z)T{) - (K4(Z)) 


The values of K,(Z), K,(Z), K3(Z) and K4(Z) are shown in table 1!.b. and the values of Ty Vv’, 
dT, and T; versus altitude, (Z), are obtained from each individual sounding. 
dt 
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oe 7+. List of Symbols and Values 


4 The following symbols have been used in the foregoing expressions. The values of some of the 
constants and also several of the mount/thermistor dimensions which have been used in the Cas! 
“. computation of data for tables 1.a and 1.b are shown. Where the values differ for the flat plate Nr 
ee and the loop mounted thermistor, the system for which the value is applicable is indicated by ‘he rege 
$ letters F and L, respectively. at 
'" wists 
| Symbol Definition Value ae 
= @ speed of sound anes 
A aerodynamic heating term es 
i Ag surface area of film 3.944 cm? meer 
A, surface area of thermistor 2.92 x 10° cm’ ed 
Alb. albedo 35 percent assumed 
> B time lag term 
Cf specific heat of film 0.2 cal gm™ (°K)” 
Cp specific heat of air at constant pressure 0.24 cal gm™ (°K)" o 
oe roa ones 
a e Cc ohmic and rf heating term ieee 
d diameter of bead 2.54 x 10° cm oe : 
D long wave emission term eae 
¥ . oN og 
ik Gir area of film multiplied by a shape factor 3.944 x 0.5 cm” S| 
aS subject to long-wave (terrestrial) radiation see 
Gi area of thermistor multiplied by a shape 2.923 x 10° x 0.5 cm’ - a 
factor subject to long-wave (terrestrial) Deets 
radiation Se 
® ah 
: Gee area of film multiplied by a shape factor 3.944 x 0.25 cm’ wage 
: supplied to short-wave (solar) radiation = 
= Gy area of thermistor multiplied by a shape 5.09 x 10° x 0.25 cm’ Bees 
as factor subject to short-wave (solar) radiation a es 
" he convective heat-transfer coefficient for film See Bibl. 4 Ae 
- hy convective heat-transfer coefficient for See Bibl. 4 
— thermistor 
Lon =z. (hplocal local convective heat-transfer coefficient for See Bibl. 4 =o 
are ‘ote film are 
re aed ay ae 
»"y o> 
b Parse 
eS. oe, 
re A-109 roe 
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Definition 

solar constant 

conductivity of lead wires 

dynamic heating correction term 

lag correction term 

radiation loss correction term 

radiation and electrical heating correction term 
Mach number 

mass of film 

ventilation velocity of sensor system, m sec” 
weighting factor 

electrical heating of thermistor due to 
measuring current and radio-frequency 
radiation 

length of lead wires 

altitude meters 


absorptivity 


absorptivity of long-wave (terrestrial) radiation 
of film 


absorptivity of long-wave (terrestrial) radiation 
of thermistor 


absorptivity of short-wave (solar) radiation of 
film 


absorptivity of short-wave (solar) radiation of 
thermistor 


cross sectional area of lead wires 
emissivity of film 

emissivity of thermistor 

angle of elevation of sun from horizontal 


viscosity 


Value Dato 


7.4 x 10°cal cm™ sec” 


(°K)" 
1.4.x 10° gm 
See Bibl. 4 


9.56 x 10°’ cal sec” 


0.33 cm 


0.1 


0.1 


0.18 


0.12 


5.07 x 10° cm’ 
0.18 

0.2 

See Bibl. 4 
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See Symbol Definition Value 
i fe) density 
ef density of film 
o Stefan-Boltzman constant 1.38 x 10°? cal cm” 
sec’ (°K)~* 
m mass of thermistor See Bibl. 4 
Ng Nusselt number See Bibl. 4 
P, Prandit number See Bibl. 4 oe 
rf recovery factor for film 0.85 : : * 
tt recovery factor for thermistor See Bibl. 4 
Rq Reynolds number based on diameter See Bibl. 4 oe 
Ry Reynolds number based on length See Bibl. 4 
R, Reynolds number based on length along the 
film See Bibl. 4 
| 5 distance along the film from stagnation point See Bibl. 4 
: S; Stanton number See Bibl. 4 
Tet effective radiation temperature of atmosphere See Bibl. 4 
Tr temperature of film --. & 
Tr temperature of thermistor, degrees K : 
FT eee temperature of atmosphere corrected, degrees K 
t time, seconds 
tr thickness of film : 
u local velocity . 
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TABLE 1 ‘oe 
Table of values of K;. Kp. K3, and Ky versus altitude for the 10-mi] bean :nerm.:s:c° using 2 loop 
mount. 
z Ky K, K2 Ky 
sec? °K I (°K) 

(m) m’ (sec) °K’ Day Night 
70 x 10° 4.65 x 10” 12.77 18.51 x 10°° 23.2 7.56 
69 4.65 10.60 14.63 18.4 6.04 
68 4.64 8.92 11.83 15.0 4.93 
67 4.64 7.70 9.80 12.4 4.13 
66 4.63 6.69 8.18 10.4 3.48 
65 4.62 5.93 6.98 8.94 3.00 
64 4.61 §.27 5.98 7,70 2.59 
63 4.59 4.72 5.17 6.69 2.26 
62 4.57 4.27 4.63 5.87 2.00 
61 4.55 3.84 3.93 §.13 1.75 
60 4.52 3.5] 3.49 4.57 1 $7 
59 4.50 3.23 3.12 4.10 14] 
38 4.48 2.95 2.78 3.66 1.27 
57 4.46 272 2.49 3.29 1.15 4 
56 4.44 2.51 2.24 2.98 1.04 . 
sé 4.43 2.342 2.05 2.73 0.956 
Sd 4.42 2.148 1.87 2.5) 0.880 
53 4.4] 2.053 1.74 2.33 0.820 
§2 4.40 1.935 1.62 2.17 0.766 
S} 4.40 1.828 1.5] 2.03 0.718 
50 4.39 1.728 1.4] 1.90 0.673 
49 3.39 1.627 1.31 i 0.629 
48 4.39 1.545 1,23 1.67 0.592 
4" 4.39 1.469 1.16 1.57 0.559 
46 4.40 1.393 1.09 1.48 0.528 
48 4.4 1.330 1.03 1.40 0.500 
44 4.4] 1.271 9.76 x 10" 1.33 0.476 
33 4.42 1.212 9.24 1.26 0.452 
42 4.42 1.160 8.79 1.20 0.430 
41 4.43 L111 8.35 LI 0.410 
40 4.43 1,065 7.95 1.09 0.392 
39 4.44 1.019 7.53 1.03 0.372 
38 4.44 0.978 7.18 0.986 0.355 
3 4.45 0.940 6.85 0.943 0.341 
36 4.45 0.903 6.52 0.897 0.323 
35 x 10° 3.46 x 107 0.868 6.22 x 10" 0.858 0.33) 
34 4.46 0.837 5.96 0.822 0.298 
33 4.47 0.804 5.68 0.785 0.285 
32 4.47 0.776 5.44 0.783 0.274 
31 4.48 0.749 5.22 0.722 0.263 
30 4.49 0.724 5.00 0.694 0.253 
29 4.49 0.700 4.79 0.666 0.243 


(m) 
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TABLE 1 (Continued) 


ky K) 
sec’°K 

m (sec) 
4.50 0.677 
4.50 0.656 
4.51 0.635 
4.52 0.616 
4.52 0.597 
4.52 0.580 
4.53 0.565 
4.54 0.550 
4.55 0.635 


4.60 
4.42 
4.25 
4.09 
3.94 
3.81 
3.69 
3.58 
3.47 


°K? 


Day 


0.640 
0.615 
0.593 
0.572 
0.551 
0.533 
0.518 
0.502 
0.487 


(°K) 


dT. 
Tatm = Tr 7 (Ky(2)V?) + kw] + (K3(Z)T*) - (Ky) 
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Night 


0.235 
0.225 
0.218 
0.210 
0.203 
0.197 
0.191 
0.185 
0.180 
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SAMPLE COPY OF ROCOB DATA ees 
(SEE ATTACHMENT 5 FOR INPUT DATA) oe 
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Sh es 
ie wy, ROCKETSONDE MESSAGE CODE 
a 
SECTION A - INTRODUCTION 
DS. 1.1 A ROCOB report or a bulletin of ROCOB reports from a land rocketsonde 
station is identified by 4M MM, = RRXX. The identifier for ROCOB SHIP 
a is SSXX. (The name ROCOB or ROCOB SHIP shall not be included in the re- 
x port.) 
1.2 The code ferm is divided into three sections as follows: 
SECTION 2 - ‘tdentification data 
SECTION 2 - Data for specified geometric altitudes 
SECTION 3 - Data for isobaric surfaces (optional) 
(Sections 2 and 3 cannot be transmitted as a separate 
report.) 
SECTION B - MESSAGE FORMAT 
a 1.1 FM39.F ROCOB — Upper-level temperature and wind (possibly air density) 
e report from land rocketsonde station. 
a SECTION 1 M,4,MLM, HCG MMJJJ = ITiii ACEC, Tee cm 
a 
iS SECTION 2 HAZ,.TT ddfff (dP PP) 
HHZ,.TT adfff (94 P,P)P)) 
. SECTION 3 112,.7,T, PyP hyhyh, daddfff 
C } eooeee eerece eoeese ates 
11Z_.T T PPhhh ddafff chee 
Tn nnonnn n n 
222,7,T, Pi Pyhyhih, dj4,f,f,F) 
( ] eeoeee eeeeve eceoe ee, 
2227 T PPhhh adfere me 
Tn nnnnn nnnnn a 
33.. eee ee ees Brgy 
44... So 00 eat ia 
= etc. 
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1.1 Symbolic form table for land stations 


M.M.M_M. 
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RRXX - ROCOB report from a land station 


Day of the month (GMT) on which the observation was taken. 


Time of launch in hours and nearest tens of minutes GMT 


u 


Month of year (01 = Jan, 12 = Dec. etc.) 


Year (980 = 1980, 981 = 1981, etc.) 


International Index Number of the observing station (II = 
Block Nbr and iii = Station Nbr, both given in H.C. Pub. 
No. 119.) 


= Reason for no report and ground equipment employed (Code 
Table 1). 


= Type of thermodynamic sensing equipment (Code Table 2). 


= Thermodynamic correction technique (Code Table 3). ae 


= Method of reducing data in temperature or wind group 
(Code Table 4). 


= Type of rocket-motor (Code Table 5). 
= Type of wind sensing equipment (Code Table 6). 
= Wind correction technique (Code Table 3). 


= Altitude, in kilometers, of the level for which data are 
reported. 


= Character of the temperature reported by TT (Code Ts -le 7). 


= Air temperature in whole degrees Celsius at the altitude 
given by HH. (For a temperature of -570, the coding is 
TT=57 and Z,=5.) When temperature is missing, two solidi 
(//) are reported for TT. 


= True direction, in tens of degrees, from which th? wind 
is blowing at the altitude given by HH (Code Table 8). 
(See note below fff.) 


= Windspeeds are preferred in meters/sec at the altitude wn 


given by HH. If windspeeds are given in knots, 50 will 
be added to YY. 
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NOTE: (1) When wind direction or speed are missing, use abnor 
solidi for ddfff as appropriate. ete 
(2) The thickness of the layer through which wind direc- an 
tion and speed are determined shall be 2 km for mee. 
both mandatory and significant levels; i.e., 1 km ones 
on each side of the altitude reported. rink 
pa. 


Indicator figure for the density group, 94 P PLP: 


Decimal point locator. The number of places to the left of 
the third significant figure. The decimal point must be so 


placed as to obtain the actual density in g/m? by pia? 3° 


NOTE: The third significant figure is always included in the 
value reported for symbol ag 


EXAMPLE: 


QO). 


tt 


Assume 120 g/m, the group is coded 90120 (i.e., a 
Assume 1.20 gym the group is coded 92120 (i.e., < = 2). 
Assume 0.281 g/m, the group 1s coded 93281 (i.e., a = 3). 
Assume 0.0788 g/m*, the group is coded 94788 (i.e., a,= 4). 


= Density in om/m* rounded to three significant figures at the 
altitude given by HH. 


= Indicator figures - Standard isobaric surface data follow; 


temperature (2,7 ene ZeT AT) is reported in the same manner 


11 
as in section 2, pressure P,P 
millibars, altitude (hjh 


potential hectometers. 


to PP is reported in whole 
4 aoe hohoh) is reported in geo- 


— 

ae 

= Indicator figures - Standard isobaric surface data follow; ae 
temperature is reported in the same manner as in Section 2, es 
pressure is reported in tenths of a millibar and altitude is enh 


reported in geopotential hectometers. 


= Indicator used for reporting pressure in hundredths of a 
millibar and altitude in geopotential hectometers. 


= Indicator used for reporting pressure in thousandths of a 
millibar and altitude in geopotential hectometers. 
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Indicator used for reporting pressure in ten-thousandths of a 
millibar and altitude in geopotential hectometers. 


Indicator used for reporting pressure in hundred-thousandths of 
a millibar and altitude in geopotential kilometers. 


Wind direction in tens of degrees (Code Table 8) at the 
reported isobaric surfaces. 


Windspeed in meters/sec at the reported isobaric surfaces. 
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ATTACHMENT 12 
TABLE OF THE LOGS OF MANDATORY PRESSURE LEVELS 
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LEVEL 
LOG 


LEVEL 
LOG 


LEVEL 
LOG 
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LOG 
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LOG 


LEVEL 
LOG 
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LOG OF MANDATORY PRESSURE LEVELS 


50 MB 
1.69897 


5 MB 
0.69897 


0.7 MB 
0.15490 


0.2 MB 
—0.699897 


0.07 MB 
—1.15490 


0.01 MB 
-2.0 


0.007 MB 
—2.15490 


0.001 MB 
—3.0 


0.0007 MB 
-3.15490 


0.0001 MB 
-4.0 


0.00007 MB 
~4,15490 


0.00001 MB 
~5.0 


30 MB 
1.47712 


3 MB 
0.47712 


0.5 MB 
~0.30102 


0.1 MB 
=f. 


0.05 MB 
—1.30102 


0.005 MB 
—2.30102 


0.0005 MB 
—3.30102 


0.00005 MB 
—4.30102 


20 MB 
1.30102 


2 MB 
0.30102 


0.4 MB 
—0.39794 


0.03 MB 
—1.52287 


0.003 MB 
~2.52287 


0.0003 MB 
—3.52287 


0.00003 MB 
~4,52287 
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10 MB 
1.0 


1 MB 
0.0 


0.3 MB 
—0.52287 


0.02 MB 
— 1.69897 


0.002 MB 
—2.69897 


0.0002 MB 
—3.69987 


0.00002 MB 
—4.69897 


7 MB 
0.84509 
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SAMPLE COPY OF HIGH-ALTITUDE METEOROLOGICAL DATA 
(SEE ATTACHMENT 5 FOR INPUT DATA) 
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SAMPLE COPY OF CARD OUTPUT 
(SEE ATTACHMENT § FOR INPUT DATA) 
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METEOROLOGICAL SOUNDING SYSTEM (MSS) 
DATA GENERAL NOVA-3/12 FORTRAN V 


ROCKETSONDE PROGRAM 
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ROCS 3 


Be 4 ABSTRACT 


The MSS rocketsonde reduction procedure is designed for use with transponder sondes only, such 
as the PWN-IOA. It consists of two programs, Pass | and Pass 2. 


This document is divided into two self-contained sections. 


Pass 1 


Pass 1 of the MSS rocketsonde reduction procedure reduces the wind and temperature data recorded 
by the MSS realtime program. Output options and surface data are entered by the operator from 
the display console. This pass produces quality control data, constant altitude data, and two disk files. 


Pass 2 


Pass 2 of the MSS rocketsonde reduction procedure extracts the HAMDATA levels from the 
intermediate scratch file and the co-rawinsonde file. This pass uses the 1,000-meter constant altitude 
file and the intermediate scratch file to produce the ROCOB coded message. 
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1.0 Module Description 
11 Main Routine 


Initializes the various control and conversion constants, reads in rocket data from diskette to scratch 
file, and controls the flow of data through the subroutines. 


1.2 Subroutine MIDINT 


Computes the midpoint of altitude layers and calculates uncorrected component winds and fall rates 
for the layers between data points. 


13 Subroutine BANDTC 

Computes ballistic motion corrections for correcting winds according to Eddy (Reference 2 in 
Bibliography). Temperature corrections are computed as described in FMH 10, Section E of Appendix 
L. 

1.4 Subroutine TWOKM 

Smooths component wind data over 2-km thicknesses, producing 2-km mean thickness winds. 


15 Subroutine SIMQ 


Uses a least squares method of smoothing, providing a best fit of the components by generating 
a second degree polynomial equation in general form. 


1.6 Subroutine PRESCOM 


Converts geometric altitude to geopotential units, computes mean virtual temperature, and computes 
atmospheric pressure using a standard form of the hypsometric equation. 


1.7 Subroutine ONEKILO 

Rearranges data in a 1-km array and computes maximum altitude recorded by the sounding. 
1.8 Subroutine OUTPUT 

Computes pressure, density, and velocity of sound, and formats tabular data for output. 

1.9 Subroutine PLOTCH 

Formats temperature data for the printer; plots temperature versus altitude as an aid to quality control. 


1.10 Subroutine COMPXYZ 


Computes X (east-west component), Y (north-south component), and Z (altitude) as well as an 
elevation angle correction for the Earth’s curvature. 


1.11 Subroutine GRAVITY 


Computes the ratio of the station’s acceleration of gravity to the acceleration of gravity at 45 degrees 
latitude, using the expression given at the top of Table 168, Smithsonian Meteorological Tables. List 
1968. 


eee 


1.12 Subroutine COMPRSR 


Computes and tabulates the differences in temperatures from the co-rawinsonde observation and 
the rocketsonde observation in the overlap region. 


1.13 Subroutine CHKOUT 


Computes and tabulates the differences in the parameters between data levels as an aid to quality 
control. 


1.14 Function IGETC 

Retrieves a character from a word or an array of words. 

1.15 Subroutine PUTC 

Inserts a character in a word or an array of words. 

1.16 Function WDIR 

Computes wind direction in meteorological polar coordinate system. 

1.17 Function WSPD 

Computes windspeed using velocity components as input data. 

1.18 Function VALUE 

Interpolates for a value or sets missing data to nines. 

2.0 MATHEMATICAL DESCRIPTION 

2.1 Main Routine 

None 

2.2. Subroutine MIDINT 

2.2.1 Computation of Mean Temperature of Two Adjacent Levels 
TP = (T + TA)/2 


where TP = mean temperature (degrees K) 


4 
it 


temperature, upper level (degrees K) 
TA = temperature, lower level (degrees K) 
2.2.2 Computation of Mean Height of Two Adjacent Levels 
H = (Z + ZA)/2 
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be 
where H = mean height (meters) 
= height, upper level (meters) 
ZA = height, lower level (meters) 


2.2.3 Computation of Mean Time of Two Adjacent Levels 
GG = (G + GA)/2 


= mean time (seconds after launch) 


= 
= 
a 
@) 
(9) 
i 


0) 
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time, upper level (seconds after launch) 
GA = time, lower level (seconds after launch) 


2.2.4 Computation of East-West and North-South Velocity Components 


(X — XA)/(GA — G) 


It 


(Y - YA)/(GA —- G) 
where Vx = east-west velocity component (meters/second) 
Vy = north-south velocity component (meters/second) 
X = east-west position component, upper level (meters) 
XA = east-west position component, lower level (meters) 
Y = north-south position component, upper level (meters) 
YA = north-south position component, lower level (meters) 
G = time, upper level (seconds after faunch) 
GA = time, lower level (seconds after lauich) 


2.2.5 Computation of Fall Rate 


Vz = (ZA — Z)/(GA — G) 
where Vz = fall rate (meters/second) 
ZA = height, lower level (meters) 


Z = height, upper level (meters) 


time, lower level (seconds after launch) 
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time, upper level (seconds after launch) 
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2.3 Subroutine BANDTC 

2.3.1 Computation of Local Acceleration of Gravity 
Gr = GRAT - 9.80616 

where Gr = local acceleration of gravity 


GRAT = ratio of local acceleration of gravity to acceleration of gravity at 45 degrees 
latitude 


9.80616 = acceleration of gravity at 45 degrees Jatitude 
2.3.2. Correction of East-West and North-South Velocity Components 
{Vz, + Vx, - ¥x,)/(G3 - Gpl} 


Vy. = Vx. - 
Xc 7 *X2 (Vz, - Vz,)/(G3 - Gp) + Gr 


° {Vz, ° UVy, - Vy,)/(G3 - Gpi} 


Vy. = Vy. -— —— 
Yo *Y2 (Vz, ~ Vz,)/G3 - Gp) + Gr 


where Vx e corrected east-west velocity component 
. V¥o = corrected north-south velocity component 
1 = upper level 
2 = intermediate level 
ve 3 = lower level 
. ; Vx = east-west velocity component 
oe Vy = north-south velocity component 
‘a Vz = vertical velocity 
G = time after launch (seconds) s 
Gr = local acceleration of gravity oe 
° 2.3.3 Computation of Ventilation Velocity mf Se 


XX = XXo = Vx, 


YY = V¥c - Vy, 


where V= 


XX,YY = 


2.3.4 Aerodynamic 


where AA = 


Ki(z) = 
V= 
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ventilation velocity 

corrected east-west velocity component 
corrected north-south velocity component 
east-west velocity component, intermediate level 


north-south velocity component, intermediate level 


= vertical velocity, intermediate level 


intermediate variables 

Heating Correction to Temperature 

Kyz) + Vv 

aerodynamic heating correction 

20—70 kilometers 

aerodynamic heating term vs altitude (see paragraph 7.3) 


ventilation velocity 


2.3.5 Time Lag Correction to Temperature 


AB = 


where AB = 


Ox 


K(z) - ((T3 rs T))/G3 te G))) 
time lag correction 


20—70 kilometers 


= time lag term vs altitude (see paragraph 7.3) 


temperature, upper level (degrees K) 
temperature, lower level, (degrees K) 
time, upper level (seconds after launch) 


time, lower level (seconds after launch) 


2.3.6 Radiation Heat Loss Correction to Temperature 


AD = 


where AD = 


Z= 


, 4 
Kyz) ° 12 
radiation heat loss correction 


20—70 kilometers 


—~ 


eo. 


K3(Z) = 
T2 = 


radiation heat loss term vs altitude (see paragraph 7.3) 


temperature, intermediate level (degrees K) 


2.3.7 Radiation and Electrical Correction to Temperature 


AC = 


Kaz) = 


Kaz) 
radiation and electrical correction 


20—70 kilometers 


radiation and electrical heating term vs altitude (variable between day and 
night; see paragraph 7.3) 


2.3.8 Computation of Corrected Temperature 


T 


where T= 


AD = 
2.4 Subroutine 
None 
2.5 Subroutine 
None 


2.6 Subroutine 


T, - AA + AB- AC + AD 

corrected temperature (degrees K) 
temperature, intermediate level (degrees K) 
aerodynamic heating correction 

time lag correction 

radiation and electrical correction 


radiation heat loss correction 


TWOKM 


SIMQ 


PRESCOM 


2.6.1 Computation of Geopotential Height 


where H = 


GRAT = 


GRAT - [(Re - Z)/(Re + Z)] 
geopotential height (meters) 


ratio of local acceleration of gravity to acceleration of gravity at 45 degrees 
latitude 


mean radius of the Earth (meters) 


geometric height (meters) 
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Ss 2.6.2 Computation of Mean Virtual Temperature ce ee 
a TV = (T + Tp)/2 noel 
a : ATA 

where TV = mean virtual temperature (degrees K) aan 

tM 

T = temperature of the level (degrees K) cat f 

| Tp = temperature, previous level (degrees K) ae 


2.6.3 Computation of the Log of the Barometric Pressure 

2 P = Pp - ((H — Hp)/(PC - TV] e 
i where P = log of the barometric pressure ee 
Pp = log of the barometric pressure, previous level 


H = geopotential height of the level 


e Hp = geopotential height, previous level Seen 
2 PC = 67.442, a constant in the hypsometric equation when units are meters ee 
TV = mean virtual temperature 

‘a 2.7 Subroutine ONEKILO a, se 

bs 2.7.1 Computation of Initial Output Height : ae 
’ HS = ((H1/1000) - 1000) + Hyj ron 
i where HS = initial output height (meters or feet) ae s 
HI = lowest input height (meters or feet) = . 
z Hint = Output interval Date 
e 2.7.2. Computation of Interpolation Ratio a ® 
: K = (HS ~ HB)/(HA ~ HB) = 

os where K = interpolation ratio 2 

- HS = output level (meters or feet) - 

s HA = upper level height (meters or feet) 


HB = lower level height (meters or feet) 
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2.7.3 Interpolation for East-West Velocity Component, North-South Velocity Component, 
Temperature, Uncorrected Temperature, and Fall Rate 


X = K - (XA — XB) + XB 


a where X = parameter interpolated 


w 
Il 


interpolation ratio 


5 


lower level parameter 


XB 


upper level parameter 
2.7.4 Interpolation of the Log of Pressure 


P = K - (PA — PB) + PB 


where P = interpolated log of pressure 
~~ K = interpolation ratio 
PA = log of pressure, lower level 
: PB = log of pressure, upper level 
: 2.8 Subroutine OUTPUT 
ae 
- 2.8.1 Computation of Barometric Pressure 
a p = oPL 
‘S 
where P = barometric pressure (millibars) 
Z PL = log of barometric pressure 
z 2.8.3 Computation of Density 
Pe @ = 348.38 - (P/T) 
oe . 
° where @ = density, grams per cubic meter 
P = barometric pressure, millibars — 
T = temperature, degrees Kelvin 
e 


348.38 = gas constant for dry air and conversion factors, with pressure in millibars 
and density in grams/cubic meter 


2.8.3 Computation of Velocity of Sound 


My VS = 331.45 - \/T/273.15 
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where VS 


velocity of sound (meters/second) 


T = temperature (degrees K) 


273.15 = freezing point of water (degrees K) 
331.45 


velocity of sound at 273.15°K (meters/second) 


2.8.4 Computation of Vertical Wind Shear 


V (Vx > Vxp) +Vy = Vyp) 


HS — HSP 


= 
”n 
H 


where WS = vertical wind shear (/second) 
Vx = corrected east-west velocity component (meters/second) 
Vy = corrected north-south velocity component (meters/second) 
VXp = corrected east-west velocity component, previous level (meters/second) 
V¥p = corrected north-south velocity component, previous level (meters/second) 
HS = output height (meters) 
HSp = previous output height (meters) 
2.9 Subroutine PLOTCH 
None 
2.10 Subroutine COMPXYZ 
2.10.1 Conversion of Range to Meters 
R = RR - 0.9144 
where R = range (meters) 
RR = range (yards) 
0.9144 = meters in a yard 


2.10.2 Computation of Geometric Height 


z {((Re + HA)’ + R’ + 2R - (Re + HA) - sin 4] — Re 


where Z = geometric height (meters) 
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a < 
oe Re = mean radius of the Earth (meters) e 

HA = station height (meters) ieee 

. R = slant range (meters) pases 
S : cones 
-" % = elevation angle pen 

2.10.3 Correction of Elevation Angle for Earth’s Curvature a =f 
¢c = 4 + {IZ - (cos $/sin $)}/2.2 - Re} c 
where $c = elevation angle corrected for the Earth’s curvature 
Fy $ = elevation angle, uncorrected 
Z = height (meters) 
Re = mean radius of the Earth (meters) 
e 2.10.4 Computation of East-West and North-South Position Components : 
7 X = R - cos ($c) - sin (8) " 
Y = R = cos ($c) - cos (6) : 
e where X = east-west position component ast 
. Y = north-south position component 
R = slant range (meters) 
"; $c = elevation angle corrected for the Earth’s curvature a 
| 6 = azimuth angle - 
2.11 Subroutine GRAVITY See 
2.11.1 Computation of Ratio of Local Gravity to Gravity at 45° Latitude ot 
GRAT = | — [00026373 - cos (2 - LAT)] + {0.000059 - [cos (2 - LAT)}*} ee 
where GRAT = ratio of local gravity to gravity at 45 degrees Jatirude 
. LAT = station latitude (degrees) 


0.0026373 constants, Smithsonian Meteorology Tables, page 488, 
0.0000059 equation 1, List 195] 


wr. toad 
ROCS 3 oe 
~ 2.12 Subroutine COMPRSR Se 
a 2.12.1 Computation of Interpolation Ratio me 
E K = (ZA ~ Z)/(ZA ~ ZB) ae 
: where K = interpolation ratio oat 
a Z = desired output height ad 
. ZA = height of the upper level a 
ZB = height of the lower level pons 
2.12.2 Interpolation for Temperature 
T = K - (TA — TB) + TB 
a where T = comparison temperature 
eo. K = interpolation ratio 


TA 


upper level temperature 


Ul 


TB lower level temperature 
2.12.3 Computation of Temperature Difference 


DIFF = TRO - TRA 


ed where OIFF = temperature difference (degrees C) 

TRO = rocketsonde temperature (degrees C) = 
cS TRA = rawinsonde temperature (degrees C) 2 
; 2.13 Subroutine CHKOUT ne 
re 2.13.1 Computation of Time Difference Between Levels in Fractional Parts of a Minute oars 
= AG = 60/(G - Gl) oe 
a where AG = time difference a 
a G = time, lower level (seconds after launch) >a 


Gl = time, upper level (seconds after launch) 
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2.13.2 Computation of Azimuth, Elevation, Range, and Height Differences for Each Minute 


ave 
yy “s 


A@ = (6 - 61) - AG 
Ae = ( - $1) + AG 
AR = (R- RI) - AG 


i} 
‘ se oe, 7 
——" 


oa eae i 


AZ = (Z - Zl) + AG 


a 


q 
a 
Je 


where A@ = azimuth difference / minute oo : 


woe 


@ = azimuth, lower level 
61 = azimuth, upper level 
AG = time difference between levels 
A¢ = elevation difference /minute 
> = elevation, lower level 
$1 = elevation, upper level 
AR = range difference / minute 
R = range, lower level 
Rl = range, upper level 
AZ = height difference / minute 
Z = height, lower level 


Z\ = height, upper level 


2.14 Function IGETC 


? 
' 
‘ 
fee 
Lite 


None 
2.15 Subroutine PUTC 


None 


fee eer ast Sait as Bie eet an 


2.16 Function WDIR 

2.16.1 Computation of Wind Direction < 
D = tan’ (Y/X) - 57.29578 

If xX < 0, 

then WDIR = 270 - D 

If xX > 0, 


If X = Oand Y <0 

. then WDIR = 180 
‘ where WDIR = wind direction, meteorological polar coordinates 
: D = wind direction, Cartesian coordinates 

X = east-west velocity vector 

Y = north-south velocity vector a 

57.29578 = degrees in a radian 
2.17 Function WSPD 
2.17.1 Computation of Windspeed 
wWsPD = x? + Y 

where WSPD = windspeed 


X = east-west velocity vector 


Y 


north-south velocity vector 

2.18 Function VALUE 

None 

3.0 INPUT 

Raw data input is from the floppy disk which has been recorded and smoothed by the realtime 


° system. Options are input in conversational mode from the display console. The co-rawinsonde 
a data input is from the floppy disk output of the rawinsonde reduction program. 


: we 
2 yoke 
$ ra fe 
Pe. res 
rs fae 
a 3 Header Record eros 
A ae eee 
‘: Words Contents : 
= 1—5 Flight identification (10 characters) 
6 Encoded sonde type 
oS 7 Number of points smoothing 
A 8 Smoothed data output rate 
tg Hard-copy option 
" 10 Hard-copy output rate 
r : U1 Scrolled data output rate 
- 12—13 Station altitude (feet) te 
. 14—15 Station temperature (degrees C) ‘ : 
‘* 16—17 Station pressure (millibars) 
18—20 Sonde identification (6 characters) 
r 21—23 Channel 1 identification (6 characters) aca 
24—25 Calibration constant 1 for channel 1 po es 
26—27 Calibration constant 2 for channel 1 setae 
28—29 Calibration constant 3 for channel | * 
30—-32 Channel 3 identification (6 characters) 
33—34 Calibration constant 1 for channel 3 mets 
e 35—36 Calibration constant 2 for channel 3 eae 
" 37—38 Calibration constant 3 for channel 3 Raeee 
- 39—40 Calibration constant for channel 2 RR 
f 41 Hygristor selection Rees 
ee 42—89 Temperature log calculation table vase 
90—137 Temperature ratio calculation table i: 
138—140 Q9 sonde identification (6 characters) 
5 : 141—142 Q9 baseline temperature ordinate vied 
oh 143144 —_Q9 baseline relative humidity Mees 
Ms 145—146  Q9 baseline temperature (degrees C) . 
wT 147 Q9 hygristor selection oe 
2 148 Flight date ne 
- 149—152 Flight starting time (hours, minutes, seconds, tenths of seconds) ae “ 
. 153 Number of blocks written eu 
oe. 154—256 Unused (zeros) 
a 
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Data Record 


Data records are packed with eight 32-word records into one 256-word block. 


Word 


Contents 


Software status word 

Time (tenths of seconds) 
Hardware status word | 

Hardware status word 2 

X position of instrument (yards) 
Y position of instrument (yards) 

Z position of instrument (yards) 
Slant range (yards) 

Ground range (yards) 

Altitude (yards above MSL) 
Azimuth (degrees from true north) 
Elevation (degrees) 

Temperature channel | (degrees C) 
Temperature channel 2 (degrees C) 
Temperature channel 3 (degrees C) 
Unused (zeros) 


Sample Options Dialogue 


Prompt 


WAS THIS A NIGHT LAUNCH (YES/NO) 


OUTPUT TO CONSOLE OR PRINTER (CNS/PRT) 


ENTER STATION ID FROM M.D.G. 


WHAT IS THE FIRST GOOD TIME 


Word 


two =— 


> Ww 


Co-rawinsonde Data Format 
Control File 
Contents 


Units switch (feet / meters) 
Winds in feet/second switch 
Interpolation switch 


Output device indicator 


’ 
Operator Input 
NO 
PRT 
01 
192 
Pate 


Control File (Continued) 
Contents 


Optical index of refraction units switch 
Wind shear in knots switch 

Surface wind x component 

Surface wind y component 

Surface wind direction 

Surface windspeed 

Station number 

Day 

Month 

Year 

Altitude increment (meters) 

Release time 

Tropopause wind direction (degrees from true north) 
Tropopause time (seconds) 

Tropopause geopotential height (meters) 
Tropopause dewpoint (degrees C) 
Tropopause geometric height (meters) 
Tropopause refractive index (optical /N) 
Tropopause pressure (millibars) 
Tropopause temperature (degrees C) 
Tropopause relative humidity (percent) 
Tropopause windspeed (meters/second) 
Station latitude (double precision) (degrees) 
Termination height (geometric meters) 
Termination pressure (millibars) 

Flight identification 

Sonde type 


Significant File 
Contents 


Altitude (meters, geometric) 
Wind direction (degrees) 
Windspeed (meters/second) one 
Temperature (degrees C) wee 
Dewpoint (degrees C) 
Pressure (millibars) 
Refractive index (N units) 
Relative humidity (percent) 


a A$ 
. 


Words 


1—2 
3—4 
5—6 
7—8 
9—10 
11—12 
13—14 
15—16 
17—18 
19—20 
21—22 
23-24 
25—26 
27—28 


40 OUTPUT 


Time 
Azimuth 
Elevation 
Slant Range 
Height 
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Mandatory File 
Contents 


Altitude (meters, geopotential) 
Wind direction (degrees) 
Windspeed (meters/second) 
Temperature (degrees C) 
Dewpoint (degrees C) 
Pressure (millibars) 

Relative humidity (percent) 


Tabulation File 
Contents 


Altitude (meters, geometric) 

Wind direction (degrees) 
Windspeed (meters/second) 
Temperature (degrees C) 
Dewpoint (degrees C) 

Pressure (millibars) 

Relative humidity (percent) 
Absolute humidity ( grams/cubic meter ) 
Density ( grams/cubic meter ) 
Index of refraction (N) 

Velocity of sound (meters/second) 
Shear velocity (seconds) 

Vapor pressure (millibars) 
Precipitable water (millimeters ) 


The output is in two forms. One is quality control outputs to the printer for use by the operator. 
The second is constant altitude data in kilometers, 250-meter intervals, and thousands of feet. 


The quality control outputs are as follows: 


Data Comparison by Level 


Seconds after midnight 
Degrees from true north 
Degrees from horizontal 
Meters from tracker 


Meters above ground level 


~~ 
oe 
Ges 

( 
yo 


Data Comparison by Level (Continued) 


Hour 

Minute 

Seconds 

Time Difference 
Azimuth Difference 
Elevation Difference 
Slant Range Difference 
Altitude Difference 


Hour of the data level (local) 
Minute of the data level (local) 
Second of the data level (local) 
Seconds between data levels 
Degrees / minute between levels 
Degrees/ minute between levels 
Meters/ minute between levels 
Meters/ minute between levels 


Temperature Comparison — Rocketsonde Vs Rawinsonde 


Altitude 

Rocketsonde Temperature 
Rawinsonde Temperature 
Temperature Difference 


The scratch file is as follows: 


Meters 

Degrees Celsius 
Degrees Celsius 
Desrees Celsius 


Words Contents 
|—2 Altitude, meters 
3—4 East-west velocity component, meters/second 
5—6 North-south velocity component, meters/second 
7—8 Uncorrected east-west velocity component, meters/second 
9—10 Uncorrected north-south velocity component, meters/second 
11—12 Temperature (degrees K) 
13—14 Uncorrected temperature (degrees K) 
1S—16 Barometric pressure (millibars) 
17—18 Fall rate, meters/second 
Constant Altitude Data 
Height Geometric feet or meters at standard intervals of | km, 250 


Wind Direction 
Windspeed 
Temperature 
Temperature Correction 
Atmospheric Pressure 
Atmospheric Density 
Fall Rate 

Speed ot Sound 

Wind Shear 


meters, and 1,000 feet 

In degrees from true north. Missing data is shown as 999 
In knots or meters/second; calm is 000 and missing is 999 
In degrees Kelvin 

In degrees Kelvin 

Millibars 

Grams/cubic meter 

In feet/second or meters/second (determined by height units) 
In knots of meters/second (determined by height units) 
Seconds" 


Ria MOS Ea tg a a eh es 


ae 
SS ROCS 3 oe 
. : NOG 
Sample Options Dialogue ie 
= Prompt Operator Input 
an Output to console or printer (CNS/PRT) PRT 
ot Plot of Temperature Vs Altitude 
) Altitude Meters ee 
Temperature Degrees Kelvin ines 
Constant Altitude Data as 
Height Geometric feet or meters at standard intervals of 1 km, 250 oi 
meters, and 1,000 feet eure 
Wind Direction In degrees from true north. Missing data is shown as 999 oe 
; Windspeed In knots or meters/second; calm is 000 and missing is 999 ae 
eo. Temperature In degrees Kelvin re 
. Temperature Correction In degrees Kelvin ae 
Atmospheric Pressure Milubars oo 
Atmospheric Density Grams/cubic meter : a rie. 
f Fall Rate In feet/second or meters/second (determined by height units) . cok 
i Speed of Sound In knots or meters/second (determined by height units) ie 
Wind Shear Seconds" oo 
A scratch file is output to diskette for Pass 2 of the program. This data is as follows: a7 
Words Contents 
m 1—2 Altitude (meters) 
: 3—4 East-west velocity component (meters/second) 
5—6 North-south velocity component (meters/second) oe 
7—8 Uncorrected east-west velocity component (meters/second) ; 
9—10 Uncorrected north-south velocity component (meters/second) 
1i—12 Temperature (degrees K) 
e 13—14 Uncorrected temperature (degrees K) . 
- 15—16 Barometric pressure (millibars) 
- 17—18 Fall rate (meters/second) 


D-AL63 634 ane STANDARDS FOR RANGE METEOROLOGICAL DATA REDUCTION 3/3 
one 2 bi agg gh RANGE slag cee WHITE 
INDS MISSILE RANGE NM_INTER-R. . 
UNCLASSIFIED TRI G-STANDARD- 3352-85-57 -= F/G 4/2 


—_—<———$—_* 


9 ye ty 


ee ee 


Mek 
1.6 
E——4 


STENSARDS. 1963 A 


=| 


——} 


ow 


oN Ob ces a 
DD d asasssi 


= il 


MICROCOPY RESOLUTION TEST CHART 
NATION AY 


es te 


3 


Ny 


ee nae | 


. 


y 


a 
eee hts? 


eo. S. 


ROCS 3 


5.0 ALGORITHMS 
5.1 Subroutine TWOKM 


Set up values that are to be used when calling SIMQ by multiplying elements of the matrix and 
assigning to proper place in symmetric coefficient matrix [A]. 


PX gy OE OER OSD 


Set up values that are to be used when calling SIMQ by multiplying elements of tae matrix and 
assigning to proper place in constant vector matrix [B]- 


By, = YB, + YCK s YK 


Using the coefficients returned from SIMQ, calculate a smoothed X and Y ere 


Sy = a XB, + aXB, + XB, 


Sy = a'YB, + aYB, + YB, 
where a=5 ne 
XB = coefficients returned from SIMQ ra 
Sx = smoothed east-west component « 
Sy = smoothed north-south component ne 


§.2 Subroutine SIMQ 


The least squares method provides a best fit of the components by generating a second degree 
polynomial equation, of general form 


y = kix? + kx + ko 


and computing an offset from the center of the 9- point range. Two matrixes are created, designated 
[A] and [B]. Symmetric coefficient matnx [A] is 


N 3X; 


ss . 2 
[A] =}2X; EX) 


=X? EX? 
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* 


and constant vector matrix [B] is 
2Y; 

{B] =|=X; Y; 

=X} Yj 
Note: <= indicates summation from | to N; N = 9. 


Basic matrix equation [A] [K] = [B] is solved by finding the inverse of [A], [A]"', and performing 
multiplications: 


[A]* [A] ([K] = [A]™ [B] 

(1] [K] = [A]” [BI 
[A] ' is calculated by SIMQ using the Gauss-Jordan inplace matrix inversion method. SIMQ returns 
the appropriate coefficients k,, k,, and ko, which the calling routine incorporates in calculating the 
appropriate values. 
6.0 FLOWCHARTS 


See the following pages. 


‘ 
7 
; 
’ 
F 
, 
. 
« 
‘ 
e, 
bd. 
~~ 
2. a af: 
." 
_ *, Ps 
* 
od 


nw 

8 

” 
eB 


6.1.1 etain Routine - Pass ] 


% 
we. 


fe 


* mat a’ 
ae 
“~> a . 
~ 
‘ 
AS 
~. OPEN 
a INITIALIZE STATION 
*. AND OPEN 10 FILE 
DATA FILES 
; READ READ 
TEMPERATURE STATION 
CORRECTION 1D 
‘ COEFFICIENTS 
a GET NIGHT cee 
‘oe LAUNCH . : 
: INFORMATION — 
os FROM OPERATOR wets 
- ae 
. : CALL te 
ad . GET LIST Fh. arp ’ 
r RDBLK 
f e DEVICE 


Am! 


FROM 


OPERATOR HEADER 


BLOCK 


GET STATION 
NUMBER 
FROM 
OPERATOR 


INITIALIZE 
BLOCK COUNT 
START TIME 
STATION ALTITUDE 


GET TYPE 
OF RUN 
FROM 

OPERATOR 


INITIALIZE 
OUTPUT 
TITLE DATA 


rere 
tee 
Pee 


if 


id 


oe Nae, 
oy 


INCREMENT 
BLOCK 
COUNTER 


UNPACK 
BUFFER 
RECORDS 


CONVERT 
TIME TO 
FLOATING 
POINT 


CONVERT 
TIME TO 
SECONOS PAST 
MIONIGHT 


6.1.2 Main Routine — Pass 1 
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YES 


CONVERT 


TEMPERATURE 
TO KELVIN 


NO 


TEMPERATURE 
Goon 
? 


RECORD 
GOOD 
? 


WRITE 
RECORD TO 
INTERMEDIATE 
FILE 


RT ST eS A 8 


SET 


TEMPERATURE 


TO MISSING 


a 
Ds 


ON, 
5 


Ba 
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4 


2 


. 
’ 

, 

rie ey 

Ge ey 
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6.1.3 Main Routine — Pass 1 


OPEN 
CO—RAWINSONDE 
CONTROL FILE 


READ 
CO-—RAWINSONDE 
CONTROL FILE 


TYPE 
ERROR 
MESSAGE 


CLOSE 
CO--RAWINSONDE 
ID FILE 


CALL 
PRESCOM 
COMPUTE 

PRESSURE 


COMPUTE 
GRAVITY 
RATIO 


CALL 
CHKOUT 


DATA 
CHECKOUT 
RUN 


MIO INTERVAL 
CALCULATION 


CALL 
BANDTC 


BALLISTIC 
MOTION 
CORRECTION 


CALL 
COMPRSR 
TEMPERATURE 
CORRECTION 


Two KM 
WIND 
SMOOTHER 


~ 


ROCS 3 


et ee ee ee ee ee ee ee 


‘@ 


6.1.4 Main Routine — Pass 1 Kant 


OUTPUT 
TABULATED 
DATA 


PLOT 
TEMPERATURE 
VS ALTITUDE 


at. e ON 

en . On 
awa 

Saf es 


bon a SE 
yk ere 
OPEN es 
PASS TWO eee 
SAVE FILE aD 
r. - 
WRITE 
DATA TO A 
SAVE FILE Ss 
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Subroutine MIDINT 


START 


READ 
DATA 
RECORD 


ASSIGN 
MAXIMUM 
ALTITUDE 


READ 
DATA 
RECORD 


YES 


NO 


COMPUTE 
MIOINTERVAL 
DATA 


WRITE A 
MID-INTERVAL 
RECORD 


DECREMENT 
DATA 


RECORD 
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6.3 Subroutine BANDTC 


Ron 
Nee? 
DEFINE 
START UNCORRECTED ie 
VARIABLES Paar 
eRe 
pg 
Ate 
ae 
READ tas 
UPPER red 
LEVEL COMPUTE 
RECORD BALLISTIC 
MOTION 
CORRECTIONS 
DEFINE 
UNCORRECTED were 
VARIABLES COMPUTE Moe 
TEMPERATURE oe 
CORRECTIONS . 
READ —o 
MIDDLE CLOSE 
LEVEL SCRATCH WRITE 
RECORD FILES CORRECTED AND 
UNCORRECTED 
DATA RECORD 
‘ 


DEFINE 


UNCORRECTED 
VARIABLES 


DECREMENT 
LEVELS OF 
INPUT DATA 


COMPUTE 
STATION 
GRAVITY 


RETURN 


soe 
eT eta ae 
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6.4 Subroutine TWOKM 


READ 
DATA TO 
ARRAYS 


WRITE 
DATA TO 
SCRATCH 

FILE 


READ 
DATA TO 
ARRAYS 


A 
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INITIALIZE 
SMOOTHING 
ARRAYS 


EAST - WEST 
COMPONENT 
DATA 


NORTH — SOUTH 


COMPONENT 
DATA 


COMPUTE 
Sv GOTHED 
DATA 


‘al 
YY 


pete ata 
a YPITE 
P SMOD7 PSU DATA 
/ TO STRATCH 
] FILE 
ane ene 
i 
Vv 
sieeshie tS 
DE CREMENT 
DATA 
ARRAYS 
Pesce sae ae 


| 
i 
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_Y 


pe 


REVERSE FILE 2 


DATA ANO 
WRITE TO 
FILE 1 


CLOSE 
SCRATCH 
FILES 


OPEN 
SCRATCH 
FILES 


RETUTN 


6.5 Subroutine SIMQ 


START 


SEARCH FOR 
MAXIMUM 
COEFFICIENT 
IN COLUMN 


TEST FOR 
PIVOT LESS 
THAN 
TOLERANCE 


INTERCHANGE 
ROWS IF 
NECESSARY 


DIVIDE 
EQUATION 
BY LEADING 
COEFFICIENT 


ELIMINATE 
VARIABLE 


RETURN 
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6.6.1 Subroutine PRFSCOM 


START ee 
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SET FALL 
RATE 
MISSING 


REAO (1) 
DATA 
PRESSURE 


BASELINE 
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6.6.2 Subroutine PRESCOM 


READ (1) 
DATA 
RECCRO 


YES 


OPEN 
SET LOG SCRATCH 
OF PRESSURE FILE 2 


MISSING 


WRITE (2) 
DATA TO 
FILE 


OPEN 
ROCKETOUT 
FILE 


COMPUTE 
LOG OF 
PRESSURE 


RETURN 


COMPUTE 
GEOPOTENTIAL 
HEIGHT 


YES 


TEMPERATURE 


MISSING 
? 


PREVIOUS 
TEMPERATURE 
MISSING 
? 


a ar 
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NO 


COMPUTE 
VIRTUAL 
TEMPERATURE 
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Subroutine GNERILO 


ie ie on 
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ao 
* 
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| SET OUTPUT 
i INTEDVAL 
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aa 


ANT 


FW 


SET 
CONSTANTS 


COMPUTE 
OUTPUT 
ALTITUDE 
(HS) 


HS 
. INPUT 
ALTITUDE 


(HA) 
? 


INCREMENT 
HS BY HINT 
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oe 
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6.7.2 Subroutine ONEKILO 


tte 
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ar od 


HS 
= INPUT 
ALT (MB) 


OUTPUT A 
DATA RECORD 


INCREMENT 
HS BY HINT 


PRECESS 
INPUT 
DATA 


COMPUTE 
INTERPOLATION 
RATIO 


INTERPOLATE 
FOR OUTPUT 
VALUES 


OUTPUT A 
DATA RECORD 


INCREMENT 


HS BY HINT 


OUTPUT A 
DATA RECORD 


COMPUTE 
MAXIMUM 
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WRITE MAX 
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RESET 
HINT AND 
CONVERSION 
FACTORS 


f ? 
| 
AY 


7 
tte 
os 
eats 
wie 
ae 


yay 
. 
? 


E 
« 
ed 
- 
. 


“eT 
. 


_ 


“? 


wand 
¥ 
"A "s 


6.8 Subroutine OUTPUT 


COMPUTE 
WIND 
Baba Mer 


WD 
SFEED 
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AND Go "esd- 


COMPUTE 
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ALT 
20600 


METERS 
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WRITE 
DATA TO 
ROCKETOUT 


DATA TO 
PRINTER 


ASSIGN 
PREVIOUS 
DATA LEVEL 
VARIABLES 


RETURN 
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wa 6.9 Subroutine PLOTCH 
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~~! 
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me 
we WRITE TITLE 
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: HEADINGS 
my TO PRINTER 
: POSITION INPUT 
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‘ LAST DATA 
= (TEM 
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oe 
ce COMPUTE 
a FIRST 
, OUTPUT 
a ALTITUDE 
. INITIALIZE 
OUTPUT 
ARRAY 


NO 
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H> 70000 
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COMPUTE 
INTERPOLATION 
RATIO 


INTERPOLATE 
FOR OUTPUT 
TEMPERATURE 


INSERT 
TEMPERATURE 
CHARACTER 
IN OUTPUT 
ARRAY 


WRITE 
OUTPUT 
ARRAY TO 
PRINTER 


DECREMENT 
OUTPUT 


ALTITUDE 


ENO OF DATA 
? 
YES 


REWIND 
INPUT 
FILE 


Yer et 
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k 
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ae 2A Pee BALE ‘- 


COMPUTE 
HEIGHT 
CORRECT 
ELEVATION 
ANGLE FGR 
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CALCULATE 
X ANDY 
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6.10 Subroutine COMPXYZ 
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6.12.1 Subroutine COMPRSR 


START 


READ 
RAWINSONDE 
HEADER 
INFORMATION 


INITIALIZE 
HEIGHT 
ARRAY 


INITIALIZE 
COUNTERS 
AND FLAGS 


POSITION INPUT 
e ROCKET FILE 
TO LAST 
DATA ITEM 


READ 
ROCKET 
DATA 
(Z AND T) 


BACKSPACE 
ROCKET 
FILE 


READ 
ROCKET 
DATA 
(21 AND T1) 


COMPUTE 
AN OUTPUT 
HEIGHT 
{HS) 


YES 


INCREMENT 
HS BY 
250 


SET 
BASELINE 
HEIGHT 
(HBL) 


WRITE 
TITLE AND 
COLUMN 
HEADINGS 


OPEN 
RAWINSONDE 
DATA 
FILE 


READ 
RAWINSONDE 
DATA 
HR, RT, RP 
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ne READ 
a RAWINSONDE 
] DATA 

u (HR1, RT1, RPI) 


6.12.2 Subroutine COMPRSR 


YES 


ROCS 2 


SET 
2=21 
T=T!1 


BACKSPACE 
ROCKET 
FILE 
TWICE 


READ 
ROCKET 
DATA 
(21 ANO T1) 


DECREMENT 
LEVEL 
COUNTER 


YES 


MISSING 
INPUT 


DATA 
? 


ves 


NO 


COMPUTE 
INTERPOLATION 
RATIO 


COMPUTE 
INTERPOLATED 
TEMPERATURE 


COMPUTE 
OUTPUT 
VALUES 


WRITE 
OUTPUT 
TO PRINTER 


hom 
or] 


. vo 
se se 
‘ 


wove 


6.12.3 Subroutine COMPRSR 


SET OUTPUT 
DATA 
MISSING 


HS < 20000 
? 


HS 


INCREMENT 
ARRAY 
COUNTER 


SET VARIABLES 
IN BASELINE 
SELECTION 

ARRAY 


INCREMENT 
HS BY 
250 
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REWIND 
ROCKET 
DATA 
FILE 


CLOSE 
RAWINSONDE 
DATA 
FILE 


INCREMENT 
JJ 


END 
QF DATA 


? 


SET 
BASELINE 
DATA 
MISSING 


ASSIGN 
BASELINE 
DATA 


RETURN 
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6.13 Subroutine CHKOUT 


WRITE TITLE 
AND COLUMN 
HEADINGS 
TO PRINTER 


INITIALIZE 
OUTPUT 
CONSTANTS 


YES 


a BOSE rare Nae tie geet Renee S| “SUN 
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CLOSE 
SCRATCH 
FILES 


RETURN 
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COMPUTE 
X,Y, AND Z 


YES FIRST 


DATA LINE 


COMPUTE 
LEVEL 
OUFFERENCES 


COMPUTE HOUR, 
MINUTE AND 
SECOND OF 
DATA LEVEL 


WRITE 
OUTPUT 
DATA TO 
PRINTER 


PRECESS 
DATA 
VALUES 
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6.14 Function IGETC 


START 


COMPUTE 
CHARACTER 
POINTER 


RETRIEVE 
CHARACTER 
FROM 
ARRAY 


RETURN 
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6.15 Subroutine PUTC 


START 


COMPUTE ieee dias 
POSITION enna 
POINTER Be 


INSERT 
CHARACTER 
IN ARRAY 


RETURN me 


ee yee ee 


6. 16 Function WDIR 


START 


SET 
WODIR = 180 


RETURN 


COMPUTE 
CARTESIAN 
WIND 
DIRECTION 


CONVERT 
TO 
METEOROLOGICAL 
POLAR DIRECTION 


RETURN 
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6.18 Function VALUE 


ROCS 3 


COMPUTE 
INTERPOLATED 
VALUE 


SET 
INTERPOLATED 
VALUE 
MISSING 


RETURN 
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7.0 


ROCS 3 


MATHEMATICAL AND PHYSICAL CONSTANTS 


Constant in hypsometric equation, Smithsonian Meteorological Tables, 


2'°, maximum number in 16 bits + 1 


Freezing point of water, degrees K 


7.1. Main Routine 
57.29578 Degrees in a radian 
67.442 
List 1968 
6371229.315 | Mean radius of the Earth, meters 
65536 
3600 Seconds in an hour 
60 Seconds in a minute 
273.15 
7.2 Subroutine MIDINT 
None 
7.3 Subroutine BANDTC 


The following table lists values of K,, Ky, K3, and Ky versus altitude for the 10-mil bead thermistor 


9.80616 


Acceleration of gravity at 45 degrees latitude, meters / second’ 


using a loop mount. 


K] 
() 
m 


4.65 x 10“ 
4.65 
4.64 
4.64 
4.63 
4.62 
4.61 
4.59 
4.57 
4.55 
4.52 
4.50 


TABLE 1 


K2 K3 
(=) 
(s) °K” 
12.77 18.51 x 10° 
10.60 14.63 
8.92 11.83 
7.70 9.80 
6.69 8.18 
5.93 6.98 
Scat 5.98 
4.72 5.17 
4.27 4.53 
3.84 3.93 
3.51 3.49 
3.23 3.12 


Seahs 
a a a 


ae 
‘ ROCS 3 
; ans 
wet 
ree, 
TABLE 1 (Continued) “oe 
a 
: Ki 2 K3 K4 = 
g Z 20 K ° Soe 
(: K ) ( ) (°K) pe 
(m) m’ (s) °K? Day Night a 
; eat 
a 58 x 10° 4.48 x 10~ 2.95 2.78 x 10°° 3.66 1.27 ies 
57 4.46 2.72 2.49 3.29 1.15 pees 
56 4.44 2.51 2.24 2.98 1.04 “ete 
55 4.43 2.342 2.05 2.73 0.956 ee 
54 4.42 2.148 1.87 2.51 0.880 
53 4.4] 2.053 1.74 2.33 0.820 poe 
ne 52 4.40 1.935 1.62 217 0.766 a 
51 4.40 1.828 1.51 2.03 0.718 ne 
50 4,39 1.728 1.41 1.90 0.673 Pane 
49 4.39 1.627 1.31 1.77 0.629 5 
48 4.39 1.545 1.23 1.67 0.592 
47 4.39 1.469 1.16 1.57 0.559 tie 
46 4.40 1,393 1.09 1.48 0.528 one's 
J 45 4.41 1.330 1.03 1.40 0.500 ee 
: 44 4.41 1.271 9.76 x 10" 1.33 0.476 oa 
43 4.42 1.212 9.24 1.26 0.452 pa 
42 4.42 1.160 8.79 1.20 0.430 a 
41 4.43 Li 8.35 1.14 0.410 oe 
40 4.43 1.065 7.95 1.09 0.392 ot 
a 39 4.44 1.019 7.53 1.03 0.372 oe 
38 4.44 0.978 7.18 0.986 0.355 one 
37 4.45 0.940 6.85 0.943 0.341 Dice 
36 4.45 0.903 6.52 0.897 0.323 wate 
35 4.46 0.868 6.22 0.858 0.311 ae 
34 4.46 0.837 5.96 0.822 0.298 ey 
33 4.47 0.804 5.68 0.785 0.285 a 
32 4.47 0.776 5.44 0.753 0.274 oo 
31 4.48 0.749 5.22 0.722 0.263 
30 4.49 0.724 5.00 0.694 0.253 
29 4.49 0.700 4,79 0.666 0.243 vor 
: 28 4.50 0.677 4.60 0.640 0.235 ais 
.. 27 4.50 0.656 4.42 0.615 0.225 mee 
26 4.51 0.635 4,25 0.593 0.218 Sag! 
2 4.52 0.616 4.09 0.572 0.210 
4 4.52 0.597 3.94 0.551 0.203 
23 4.52 0.580 3.81 0.533 0.197 ae 
22 4,53 0.565 3.69 0.518 0.191 oe 
Mae 21 4.54 0.550 3.58 0.502 0.185 no 
- 20 4.55 0.535 3.47 0.487 0.180 oat 
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ee Eg RE Te SG MRE MI ES Sig, lee See SMS te ee MS ae Oie, Rae came emg ea ae Nm Me cae tm de Oi a RE ae, 


a ROCS 3 mas 

: mes ae 

= 7.4 Subroutine TWOKM eee 

Ss 

None pone 

75 Subroutine SIMQ nae 

None Seat 

KAN 

7.6 Subroutine PRESCOM . ee 

None . : 

7.7 Subroutine ONEKILO ote 

1000 Output interval, meters and feet me 

250 Output interval, meters on 

7.8 Subroutine OUTPUT ae 

348.38 Gas constant for dry air with pressure in millibars and density in : ~ 

grams/cubic meter, combined with conversion factors ice 

331.45 Speed of sound at zero degree Celsius, meters/second £ : 2 

273.15 Freezing point of water, degrees K ‘.. ronan 

7.9 Subroutine PLOTCH coe 

: 1000 Output interval, meters aoe 

= 70000 Maximum output altitude fe 

s 7.10 Subroutine COMPXYZ rere 
% 0.9144 Meters/yard 

‘ 2.2 Constant used in correcting elevation angle for the Earth’s curvature eee 

7.11 Subroutine GRAVITY | 

57.29578 | Degrees in a radian a 

° 0.0026373 Gravity ratio constants from Smithsonian Meteorological Tables, chee 

0.0000059 page 488, equation 1, List 1951 ne 

7.12 Subroutine COMPRSR on 


0.3048 Conversion factor for meters to feet 


273.15 Freezing point of water, degrees K 
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1.0 Module Description 
1.1 Main Routine 


Initializes the various control ane conversion constants, and controls the flow of data through the 
subroutines. 


1.2 Subroutine ETFORM 


Formats data for entry into the high altitude data formatting routine to develop data for archiving 
and publication. 


13 Subroutine CODIT 

Final formatting and output routine for data developed in Subroutine ETFORM. 

1.4 Subroutine RDRAWIN 

Reads co-rawinsonde data for input into the archiving and publication formatting program. 
1.5 Subroutine ROCOB 


Encodes all rocketsonde parameters into the international ROCOB code and stores data for later 
transmission. 


1.6 Subroutine MANDL 

Mandatory data level inserter called by ROCOB to place data levels into the proper place in the code. 
1.7 Subroutine LEFTZERO 

Inserts leading zeros on integer numbers when required for output. 
1.8 Subroutine TSIGN 

Converts a decimal integer to a positive or negative character. 

1.9 Function VALUE 

Interpolates for a value or sets missing data to nines. 

1.10 Function WDIR 

Computes wind direction in meteorological polar coordinate system. 
1.11 Function WSPD 


Computes wind velocity using velocity components as input data. 
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2.0 MATHEMATICAL DESCRIPTION 


2.1 Main Routine 

None 

2.2 Subroutine ETFORM 

2.2.1 Conversion of Baseline Height to Nearest Geopotential Decameter 


H = {[GRAT - (Re - Z)/(Re + Z)] + 5}/10 


where H = baseline height, geopotential decameters 


GRAT = ratio of local acceleration of gravity to acceleration of gravity at 45 
degrees latitude 


Re = mean radius of the Earth, meters 
Z = geometric height, meters 


5 = rounding factor 


10 converts meters to decameters 


2.2.2 Computation of Baseline Pressure 


p= (10° . 100) + 0.5 
where P = baseline pressure 
PL = log of baseline pressure 
100 = formats pressure for output 
0.5 = rounding factor 


2.2.3 Computation of an Interpolation Ratio for Significant Level Checking 


K = (H — HB)/(HA — HB) 
where K = interpolation ratio 
H = height of level being checked 
HA = height, upper level 


HB = height, lower level 
2.2.4 Interpolation for Significant Level Test Temperature 
ST = K - (TA — TB) + TB 
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significant level test temperature, degrees Kelvin 
interpolation ratio “age 
temperature, upper level, degrees Kelvin 


temperature, lower level, degrees Kelvin 


2.2.5 Computation of an Interpolation Ratio for Obtaining Even Decameter Data 


K 


where K 


H1 


HA 


HB 


2.2.6 Conversion 


Hes 


where H 


Hl 


10 


(H1 — HB)/(HA — HB) 
interpolation ratio 

Output height, meters 
height, lower level, meters 


height, upper level, meters 


of Output Height to Decameters 


H1/10 
height, decameters 
height, meters a 


conversion factor 


2.2.7 Interpolation for Output Parameters 


X = 
Y 
T= 
P 
F 


UX = 


x KR KR KR KR KK 


UY 
where X = 


xXA = 


K - (XA — XB) + XB 


o, 
> 
| 


YB) + YB 

- (TA - TB) + TB 

- (PA — PB) + PB 

- (FA — FB) + FB 

- (UXA - UXB) + UXB 

- (UYA — UYB) + UYB 

output east-west position component 
lower level east-west position component 
upper level east-west position component 


output north-south position component 


B-156 


WTO yee Te ee UR SS ee 
SO wae stkcee ee et ee Pe Boe Neen er oe be a 


‘ YA = lower level north-south position component aos 
- : YB = upper level north-south position component Re 
4 T = output temperature ann 
2 TA = lower level temperature re 
% TB = upper level temperature be 
P = output log of pressure oe 
. PA = lower level log of pressure = 
- PB = upper level log of pressure oe 
F = output fall rate r : 
FA = lower level fall rate 
FB = upper level fall rate foe 
ze UX = output uncorrected east-west position component Se 
= UXA = lower level uncorrected east-west position component seat 
= e UXB = upper level uncorrected east-west position component ot 
UY = output uncorrected north-south position component an 
UYA = lower level uncorrected north-south position component mh 
= UYB = upper level uncorrected north-south position component ost 
K = interpolation ratio oa 
“ 2.2.8 Computation of an Interpolation Ratio for Obtaining the Constant Pressure Output Data ae 
eS K = (P - PB)/(PA ~ PB) os 
. where K = interpolation ratio ee 
: P = log of pressure, mandatory output levels 
a PA = log of pressure, lower level Doh 
PB = log of pressure, upper level wee 
2.2.9 Interpolation for Height, East-West Position Component, North-South Position Component, ae 
Uncorrected East-West Position Component, Uncorrected North-South Position Component, Rae 
Temperature, Uncorrected Temperature, and Fall Rate oe 
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R- Ka 
e X = K - (XA ~ XB) + XB od, Ge 
: Y= K - (YA- YB) + YB Bh gs 
T = K - (TA - TB) + TB rer 
P = K - (PA -~ PB) + PB we 
F = K - (FA — FB) + FB ih 
UX = K - (UXA — UXB) + UXB eae 
UY = K - (UYA = UYB) + UYB oa 
where X = output east-west position component = 
XA = lower level east-west position component a 
XB = upper level east-west position component ae 

Y = output north-south position component - 

YA = lower level north-south position component 
YB = upper level north-south position component 2 a 
T = output temperature ae 
TA = lower level temperature ee oe 
TB = upper level temperature <; ee 
P = output log of pressure a 
PA = lower level log of pressure i : 

PB = upper level log of pressure ae 
F = output fall rate “ 
FA = lower level fall rate BS 
FB = upper level fall rate — 
UX = output uncorrected east-west position component a 
UXA = lower level uncorrected east-west scallion component = 
UXB = upper level uncorrected east-west position component ie 

UY = output uncorrected north-south position component . 

UYA = lower level uncorrected north-south position component : 
UYB = upper level uncorrected north-south position component 
K = interpolation ratio a 
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: oo 2.3 Subroutine CODIT 
ee 2.3.1 Computation of Velocity of Sound 


ie VS = 331.45 - (1/273.15)'” 
> where VS = velocity of sound, meters per second 
= T = temperature, degrees Kelvin 
331.45 = speed of sound at zero degree Celsius, meters/second 


2.3.2 Computation of Barometric Pressure 


p = j0PL 
= where P = barometric pressure (millibars) 
ae PL = log of barometric pressure 
« 2.3.3 Computation of Density 
* @ = 348.38 - (P/T) 
where @ = density, grams per cubic meter 
, e P = barometric pressure, millibars 
T = temperature, degrees Kelvin 
348.38 = gas constant for dry air and conversion factors, with pressure in millibars 
and density in grams/cubic meter 
2.4 Subroutine RDRAWIN 
a None 
: 2.5 Subroutine ROCOB 
de None 
2.6 Subroutine MANDL 
None 
©. 2.7 Subroutine LEFTZERO 
Ds None 
2 2.8 Subroutine TSIGN 


= te, None 


2.9 ‘Function VALUE 
None 
2.10 - Function WDIR 
2.10.1 Computation of Wind Direction 
D = tan (Y/X) - 57.29578 


If X < 0, 

then WDIR = 270 - D 

If X > 0, 

then WDIR = 90-D 

If X = 0 and Y 20 

then WDIR = 360 

If X = 0 and Y <0 

then WDIR = 180 

where WDIR = wind direction, meteorological polar coordinates 


D = wind direction, Cartesian coordinates 
X = east-west velocity vector 
Y = north-south velocity vector 


57.29578 


degrees in a radian 
2.11 Function WSPD 

2.11.1 Computation of Windspeed 
Vx? + Y' 


windspeed 


WSPD 


u 


where WSPD 


Ml 


x 


east-west velocity vector 


Y 


i] 


north-south velocity vector 
3.0 INPUT 


Input to Pass 2 is from four sources. First is the scratch file output from Pass 1, second is for 
the constant altitude data output from Pass 1, third is the options that are input in conversational 
mode, and fourth is the co-rawinsonde data. 
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e . ae 
say 
ROCS 3 
eae 
: : aes 
aa The scratch file is as follows: waht 
4 e Pe ” reat o 
ae es 
Words Contents , 
ie a) 
ew 
- 1—2 Altitude, meters He 
. 3—4 East-west velocity component, meters/second eee 
ay 
~ 5—6 North-south velocity component, meters/second 
7—8 Uncorrected east~west velocity component, meters/second 
9—10 Uncorrected north-south velocity component, meters/second 
11—12 Temperature (degrees K) a 
13—14 Uncorrected temperature (degrees K) ee 
15—16 Barometric pressure, millibars mee 
17—18 Fall rate, meters/second - 
Constant Altitude Data men 
Height Geometric feet or meters at standard intervals of 1 km, 250 s 
meters, and 1,000 feet 
- Wind Direction In degrees from true north. Missing data is shown as 999 ps 
e Windspeed In knots or meters/second; calm is 000 and missing is 999 se 
eel 
Temperature In degrees Kelvin ols 
ee 
Temperature Correction In degrees Kelvin wo 
Atmospheric Pressure Millibars mrad 
oe Atmospheric Density Grams/cubic meter : 
oe Fall ‘Rate In feet/second or meters/second (determined by height units) ore 
Kd 
" Speed of Sound In knots or meters/second (determined by height units) 
Wind Shear Seconds” 
Sample Options Dialogue 
“- Prompt Operator Input “ 


Output to console or printer (CNS/PRT) PRT 


ad eah EET ET EEE ES MT PON A 
iat ta A 


Contents 


Units switch (feet /meters) 

Winds in feet/second switch 

Interpolation switch 

Output device indicator 

Optical index of refraction units switch 
Wind shear in knots switch 

Surface wind x component 

Surface wind y component 

Surface wind direction 

Surface windspeed 

Station number 

Day 

Month 

Year 

Altitude increment (feet) y 
Release time 

Tropopause wind direction (degrees from true north) 
Tropopause time (seconds) 

Tropopause geopotential height (feet) 
Tropopause dewpoint (degrees C) 
Tropopause geometric height (feet) 
Tropopause refractive index (optical/N) 
Tropopause pressure (millibars) 
Tropopause temperature (degrees C) 
Tropopause relative humidity (percent) 
Tropopause wind speed (knots /feet / second / meters / second ) 
Station latitude (double precision) (degrees) 
Termination height (geometric feet) 
Termination pressure (millibars) 

Flight identification 

Sonde type 


a Se Ve ee 


Sor Se Sea Le a Se ea oS Ne ee 


Contents 


Altitude (feet, geometric) 
Wind direction (degrees) 
Windspeed (meters/second) 
Temperature (degrees C) 
Dewpoint (degrees C) 
Pressure (millibars) 
Refractive index (N units) 
Relative humidity (percent) 


Mandatory File 
Contents 


Altitude (feet, geopotential) 
Wind direction (degrees) 
Windspeed (meters/second) 
Temperature (degrees C) 
Dewpoint (degrees C) 
Pressure (millibars) 

Relative humidity (percent) 


Tabulation File 
Contents 


Altitude (feet, geometric) 
Wind direction (degrees) 
Windspeed (meters/second) 
Temperature (degrees C) 

Dew point (degrees C) 

Pressure (millibars) 

Relative humidity (percent) 
Absolute humidity (grams/cubic meter ) 
Density (grams/cubic meter ) 
Index of refraction (N) 
Velocity of sound (knots) 
Shear velocity (seconds”') 
Vapor pressure (millibars) 
Precipitable water (millimeters) 
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Say 
Output consists of the HAMDATA file and the ROCOB file. The first part of the HAMDATA 
file is as follows: 

Words Contents 
l Height, geometric decameters for the significant and 

mandatory data levels, geopotential decameters for 

the constant pressure levels 
2—4 Wind direction (degrees from true north) 
5 Windspeed (meters/second) 
6 Uncorrected north-south velocity component 

(meters / second) 
7 Uncorrected east-west velocity component 

(meters / second) 
8 North-south velocity component (meters/second) 
9 East-west velocity component (meters/second) 
10 Fall rate (meters/second) 
1] Temperature (degrees C) 
12 Temperature (degrees C) \- 
13—15 Barometric pressure (millibars) 
16—18 Density (grams/cubic meter ) 
19 Velocity of sound (meters/second) 


The second part of the HAMDATA file consists of the data in the mandatory and significant parts 
of the co-rawinsonde data and is as follows: 


Mandatory File 


Contents 


Altitude (feet, geopotential) 
Wind direction (degrees) 
Windspeed (knots) 
Temperature (degrees C) 
Dewpoint (degrees C) 
Pressure (millibars) 
Relative humidity (percent) 


P08 Rg SE SS SO 


5—6 
7—8 
9—10 
11—12 


Contents 


Altitude (feet) 

Wind direction (degrees) 
Windspeed (meters/second) 
Temperature (degrees C) 
Dewpoint (degrees C) 
Pressure (millibars) 
Refractive index (N units) 


The ROCOB file is described in Federal Meteorological Handbook !0 as follows: 


ROCKETSONDE MESSAGE CODE 


SECTION A - INTRODUCTION 


1.1 A ROCOB report or a bulletin of ROCOB reports from a land rocketsonde 
station is identified by MMMM. = RRXX, The identifier for ROCOB SHIP 


is SSXX. (The name ROCOB or ROCOB SHIP shall not be included in the re- 


port.) 


1.2 The code form is divided into three sections as follows: 


SECTION 1 - Identification data 


SECTION 2 - Data for specified geometric altitudes 


SECTION 3 - Data for isobaric surfaces (optional) 


(Sections 2 and 3 cannot be transmitted as a separate 
report.) 


SECTION B - MESSAGE FORMAT 


1.1 FM39.F ROCOB — Upper-level temperature and wind (possibly air density) 
report from land rocketsonde station. 


SECTION 1 


SECTION 2 


array ar YYCG, Iliii By ODE DCM Few ew ower 
HHZ,,TT ddf ff ta ta iat ty 
HHZ, TT adfff (94 P P,P) 


SECTION 3 


MiMM.M; 


YY 


GG 
g 
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112,7,T, PAP ALA, ddfffF 
Y1Z.T T PPhhh ddffre 
Tnn nnannn nnnnn”a 
222,7,T, P,P jb pin, d,4,f,f,f) 
22200 Ty en nn qa tnen 
33... seeiers Sid.eiece 
44... Seisee. jj ~~ o reteetecc 
etc. 


SECTION C - DEFINITIONS 


1.1 Symbolic form table for land stations 


RRXX - ROCOB report from a land station 

Day of the month (GMT) on which the observation was taken. 
Time of launch in hours and nearest tens of minutes GMT 
Month of year (01 = Jan, 12 = Dec. etc.) 

Year (980 = 1980, 981 = 1981, etc.) 

International Index Number of the observing station (IT = 
Block Nbr and iii = Station Nbr, both given in H.O. Pub. 


No. 119.) 


Reason for no report and ground equipment employed (Code 
Table 1). 


Type of thermodynamic sensing equipment (Code Table 2). 


Thermodynamic correction technique (Code Table 3). 


Method of reducing data in temperature or wind group 
(Code Table 4). 


Type of rocket motor (Code Table 5). 


Type of wind sensing equipment (Code Table 6). 
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ROCS 3 
e,, = Wind correction technique (Code Table 3). 
HH = Altitude, in kilometers, of the level for which data are 
reported. 
Zn = Character of the temperature reported by TT (Code Table 7). 
TT = Air temperature in whole degrees Celsius at the altitude 
given by HH. (For a temperature of -570, the coding is 
TT=57 and Z_=5.) When temperature is missing, two solidi er 
(//) are reported for TT. ies oe 
dd = True direction, in tens of degrees, from which the wind ee 
| is blowing at the altitude given by HH (Code Table 8). ‘eae 
@ 
; (See note below fff.) on 
fEE = Windspeeds are preferred in meters/sec at the altitude 


given by BH. If windspeeds are given in knots, 50 will 
be added to YY. 


NOTE: (1) When wind direction or speed are missing, use 
solidi for ddfff as appropriate. 


(2) The thickness of the layer through which wind direc- 
tion and speed are determined shall be 2 km for 
both mandatory and significant levels; i.e., 1 km 
on each side of the altitude reported. 


9 = Indicator figure for the density group, shade ide es oe 


a = Decimal point locator. The number of places to the left of 
the third significant figure. The decimal point must be so 


placed as to obtain the actual density in g/m by p1°3?- 


NOTE: The third significant figure is always included in the hg 
value reported for symbol q. eee 


EXAMPLE : 

Assume 120 ayn the group is coded 90120 (i.e., oo = 0). 

Assume 1.20 g/m, the group is coded 92120 (i.e., qs = 2). 

Assume 0.281 g/m, the group is coded 93281 (i.e., os = 3). 

Assume 0.0788 ane the group is coded 94788 (i.e., os = 4). e 
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E PY Py Py = Density in gm/m> rounded to three significant figures at the wee 
j altitude given by HH. 


11 = Indicator figures - Standard isobaric surface data follow; 


temperature (2,7, T, din 2,77) is reported in the same manner 
ecece i i h 
17) PP.) is reported in whole 


| as in section 2, pressure P. 
millibars, altitude (hyhyh, Sieve hhh) is reported in geo- 
potential hectometers. 


4 22 = Indicator figures - Standard isobaric surface data follow; 
temperature is reported in the same manner as in Section 2, 
pressure is reported in tenths of a millibar and altitude 
reported in geopotential hectometers. 


2 33 = Indicator used for reporting pressure in hundredths of a 
2 millibar and altitude in geopotential hectometers. 
44 = Indicator used for reporting pressure in thousandths of a 


. millibar and altitude in geopotential hectometers. 


4a ‘ 55 = Indicator used for reporting pressure in ten-thousandths of a . 
millibar and altitude in geopotential hectometers. 


66 = Indicator used for reporting pressure in hundred-thousandths of 
a millibar and altitude in geopotential kilometers. 


i did = Wind direction in tens of degrees (Code Table 8) at the 
reported isobaric surfaces. 


fifif, = Windspeed in meters/sec at the reported isobaric surfaces. 
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~ "54 Subroutine MANDL 


Convert mandatory level pressure to whole number for output in Rocket code using 


= IPS 


= where IPS 
PL(J) 
Fl 
\. » 
ca 
s JS 


doPLO) + F1) (1045/10) 


pressure at mandatory level encoded as a whole number 


log of the pressure at the mandatory level being converted to 
a whole number : 


round off factor used to move decimal point; this factor is 
based on PL(J) 


If PL(J) < 0, Fl = 0.05 
< -1, Fl = 0.005 
< -2, Fl = 0.0005 
< -3, Fl = 0.00005 
< —4, Fl = 0.000005 


intermediate factor based on the pressure level indicator J 
which is set by the number of levels processed 


If J > 31, IS =6 
< 31, JS =5 
< 26, IS =4 
< 21, JS =3 
< 16, JS =2 
< 10, JS=1 


5.2 Subroutine TSIGN 


Convert the tenths value of temperature to a positive or negative character. 


The temperature is read as a positive or negative real number. If it is positive the tenths value 
is converted to a character number. If it is negative the tenths value is converted to a coded character. 


Numeric Value Character 


a HDun & WN —- © 
4 WDNrA hk WN — © 
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oo Numeric Value Character 


fe 8 8 we 

9 9 
-0 ! 
-1 J 
-2 K 
~3 L 
-4 M 
-5 N 
-6 O 
-7 P 
‘oa 8 Q 
-9 R 


53 Subroutine LEFTZERO 


Generate an array of single digit numbers from successive divisions of the input integer number. 
This routine is used to put leading zeros on integer numbers for output formatting. 


6.0 FLOWCHARTS 
See the following pages. u 
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6.1 Main Routine — Pass 2 


CALL 
ROCOB 
WORK UP 
ROCKET CODE 
RELEASE 


DISKETTE 
DRIVE 


DELETE 
PREVIOUS 


FILES 


INITIALIZE 
DISKETTE 
DRIVE 


STOP 


OPEN 
PASS 1 
SAVE FILE 


READ 
PASS 1 
SAVE FILE 


ASK 
OPERATOR 
FOR OUTPUT 
DEVICE 


e ao . 
oa ~ 
[tiv \ oe 
Cone 


WORK UP cs 
HAMDATA sat 
~ 


’ 6.2.1 Subroutine ETFORM 


START 


OPEN 
HAMDATA 
AND SCRATCH 
FILES 


INITIALIZE 
CONSTANTS 


READ A LEVEL 
OF INPUT 
DATA FROM 
SCRATCH FILE 


PRESSURE 


MISSING 
? 


COMPUTE 
BASELINE 
DATA 


WRITE 
BASELINE 
DATA TO 

HAMDATA FILE 
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INITIALIZE 
MAX AND 


REWINO 
SCRATCH 
FILE 


MIN FLAGS 
(IS? AND IS2) 


WRITE COMPUTE 


TITLE TO UPPER 
HAMDATA BOUNDARY 
FILE ALTITUDE (H1) 


INITIALIZE 
FLAGS AND 
DATA 
VARIABLES 


FIRST 
DATA LEVEL 
? 


INITIALIZE 
MAXIMUM 
AND MINIMUM 

TEMPERATURE 


° 
* 

14% tf 
soe e 
ae ee 


e 
’ 


UA 
vee 


he 
> 


man 


eT 


6.2.2 Subroutine ETFORM 


DETERMINE 
SIGNIFICANT 


LEVEL BOUNDARY 
PARAMETERS 


1S 
THERE A 

SIGNIFICANT 

LEVEL 
? 


ASSIGN 
SIGNIFICANT 
LEVEL DATA 
TO OUTPUT 
VARIABLES 


QUTPUTA 
DATA RECORD 


COMPUTE 
AN 
INTERPOLATION 
RATIO 


INTERPOLATE 
FORA 
MANDATORY 
OUTPUT LEVEL 
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OUTPUT A 
DATA RECORD 
BACKSPACE 


SCRATCH 
FILE 


TT I EF Co er a vA 
. 


BACKSPACE 
SCRATCH 
FILE TWICE 


READ 
LAST DATA 
LEVEL 


OUTPUTA 
DATA RECORD 


REWIND 
SCRATCH 
FILE 


y 
ry 
‘ 
cy 
t 
. 
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a - 70 MATHEMATICAL AND PHYSICAL CONSTANTS 
me 7.1. Main Routine os 
« 
oe Log of Mandatory Pressure Levels 
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ROCS 3 
ATTACHMENT 3 
CO-RAWINSONDE FILE (SIGNIFICANT AND MANDATORY) 
PASS 2 
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ROCS 3 
ATTACHMENT 4 
HAMDATA OUTPUT FILE 

PASS 2 
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